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Abstract 


I prove that the Bethe roots describing either the ground state or a certain class of 
"particle-hole" excited states of the XXZ spin-1/2 chain in any sector with magneti¬ 
sation m e [0; 1 /2] exist and form, in the inhnite volume limit, a dense distribution 
on a subinterval of R. The results holds for any value of the anisotropy A > -1. In 
fact, I establish an even stronger result, namely the existence of an all order asymp¬ 
totic expansion of the counting function associated with such roots. As a corollary, 
these results allow one to prove the existence and form of the inhnite volume limit of 
various observables attached to the model -the excitation energy, momentum, the zero 
temperature correlation functions, so as to name a few- that were argued earlier in the 
literature. 


Introduction 

The XXZ spin-1/2 chain refers to a system of interacting spins in one dimension described by the Hamiltonian 
L 

Ha - . (0.1) 

a=\ 

Ha is an operator on the Hilbert space of the model hxxz - with ha - C^. The matrices tr"', w = x,y,z 

are the Pauli matrices and cr^ stands for the operator on hxxz which acts as the Pauli matrix cr”' on ha and as 
the identity on all other spaces appearing in the tensor product defining hxxz- The Hamiltonian depends on two 
coupling constants : 7 > 0 which represents the so-called exchange interaction and A which takes into account 
the anisotropy in the coupling between the spins in the longitudinal and transverse directions. Finally, the chain 
consists of an even number of sites L £ 2N. 

The Hamiltonian Ha commutes with the total spin operator 

L 

a=l 
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The Hilbert space of the model hxxz decomposes into the direct sum 
L 

^xxz ^ ^ with h^z ^ {v e i)xxz '■ ■ v ^ (L - 2N) ■ vj , (0.3) 

N=0 

what turns Ha into a block diagonal operator relatively to the above decomposition. I denote by the restriction 
of Ha to every subspace h^z' 

The spectrum and eigenvectors of the isotropic limit A = 1 of the XXZ chain have been first studied by Bethe 
|2] in 1931 by means of the celebrated Bethe Ansatz. Then, in 1958, Orbach extended the approach to the case of 
the XXZ chain with a general coupling A. Within Bethe’s Ansatz, the eigenvectors and eigenvalues of H^^ 
are parametrised by a collection of N numbers {Tq}^, the so-called Bethe roots. The Bethe roots are constrained 
to satisfy a system of transcendental equations called the Bethe equations. For -1 < A < 1 when the anisotropy 
can be parametrised as A = cos(4') with ^ e]0; 7r[, these take the form 


r-r fsinhCTa^^Tfoj^'l / sinh(i^/2 -/la) \^ _ 
\ sinh(/lfo -Aa + iol v sinh(i^/2 -i- A^) I 


(0.4) 


When A = 1, the Bethe equations degenerate into the rational form 


|~r f dg - Tfa -t i ] / i/2 - 

1 dfo - /Iq-I- i J li/2 -I- /la / 


(-1)'^+!, 


(0.5) 


Finally, for A > 1, the anisotropy can be parametrised as A = cosh(^) with ^ e and the Bethe equations take 
the form 


pr f sin(/la -/Ij, -tip ) / sin(i^/2-/la) \^ ^ t-n^+i 

1=1 - 4a -I- iO I lsin(i^/2 + Aa)J 


(0.6) 


The question of the completeness of the Bethe Ansatz, namely whether the set of solutions to (10.41) . (10.51) or (10.61) . 
is in one-to-one correspondence with the set of eigenvectors of H^^ is tricky and remained open for quite a long 
time. A positive answer has been given, for the XXX chain, by Mukhin, Tarasov and Varchenko Il29l provided 
that one agrees to slightly change the perspective and to characterise states in terms of polynomial solutions to an 
appropriate T - Q equations. Also, completeness was established for certain generic inhomogeneous variants of 
the XXZ chain ll^ . 

Independently of completeness issues, there are numerous other questions related with the study of the equa¬ 
tions (I0.4I) - (I0.6I) . One relates to identifying the solution of (10.41) . (10.51) or (10.61) . depending on the value of A, giving 
rise to ’s ground state, viz- the eigenvector associated with the lowest eigenvalue. The answer has been obtai¬ 
ned by Yang-Yang in [361. By applying a variant of the Perron-Frobenius theorem, Yang and Yang showed that 
h|^^ admits a unique ground state. In order to identify the roots describing the ground state, it is convenient to 
rewrite the Bethe equations in their logarithmic form 


ZTT 




IttL 


■ 4fo) -I- 


a=l 


N+l 

2L 


4 

L 


- — with 4 £ 


where for -1 < A = cos(/;') < 1, ^ e]0 ;n'[ 


p(/l) = iln 


sinh(i^/2 -i- /l)\ 
sinh(i/;' 12- A)) 


0(A) 


/ sinh(i/f -I- T) \ 

I sinh(i^ -/!)/’ 


(0.7) 


( 0 . 8 ) 
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for A = 1 


p(/l) = iln 


ill + A\ 
ill- a) 


0(A) - iln 


i + A\ 


and, for A = cosh(^) > 1, ^ £ R"*", 


(0.9) 


/\ 

V(A)^§(A,ai) , 9{A)^§(A,0 with i?(d,77) = f. (0.10) 

J sin(ju + irj) sin(yu - i?j) 

0 

Yang and Yang proved that the ground state for A > -1 is obtained from a real valued solution to (10.71) corres¬ 
ponding to the specific choice of integers 7’^ = a. More precisely, Yang and Yang were analysing a reparametrised 
version of the Bethe equations in terms of ka variables, Aa - fAika) for some explicit /a. It was, in fact, the 
equation in terms of the variables ka that has been initially obtained by Orbach BOl . Yang and Yang showed that 
the transformed counterpart of (10.71) admits a unique real valued solution when -1 < A < 0. Furthermore, they 
showed that the transformed counterpart of (10.71) when (a - a admits, in fact, solutions for any A > -1 and that, 
among these, there exists one such that A ^a(A) is continuous in A > -1. Yang and Yang established that 
it is precisely this particular solution that gives rise to the Bethe roots parametrising the ground state of I 
will refer, henceforth, to this solution as the ground state Bethe roots. Yang and Yang, however, did not prove the 
uniqueness of solutions to (10.71) for A > 0 and ia - a. 

On top of the one parametrising the ground state of the logarithmic Bethe Ansatz equations admit many 
other solutions upon varying the choices of integers ia- For instance, when the choice of the ground state integers 
is slightly perturbed, as 


(a = a for a £ 11; AI \ {/Ja)i = pa with 


/2«e|[l; AI 

P^£|[-Ma;Ma + AI\[[1; A3 


( 0 . 11 ) 


with Ma an integer depending on A, then the real-valued solution to the logarithmic Bethe Ansatz equations -if 
they exist- define so-called particle-hole excited states. One can also have complex valued solutions to the Bethe 
Ansatz equations, as already observed by Bethe O. See |41 for an extensive numerical analysis thereof in the case 
of small length XXX chains. In the present paper, I will not discuss the complex valued solutions. 


From the point of view of practical applications one is mostly interested in the behaviour of observables 
attached to the model in the so-called thermodynamic limit L —> -i-oo. Such observables can be the energy or 
momentum of an eigenstate or some correlation function. In practical situations this amounts to computing either 
the limit or the first few terms in the large-L asymptotic expansion of sums of the type 

1 ^ 

^ 2 ]/( 4 ;) ( 0 . 12 ) 

a=l 

where / is some sufficiently regular function, {Aa\^ are the Bethe roots describing the ground state or some exited 
state "close" to it and the integer A labelling the spin sector to which the Bethe vector belongs is L dependent and 
grows with L in such a way that A/L ^ D £ [0; 1/2]. 

In 1371, Yang and Yang affirmed that the limit exists in the case of the Bethe roots for the ground state and 
that, for any sufficiently regular function /, it holds 


lim 

L-^+oo 


{iZw} - ff(A)p(Alq)-dA. 


(0.13) 
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The pair {q,p{A \ q)) appearing above corresponds to the unique solution to the system of equations for the 
unknowns {Q,p{X \ Q)) : 


Q 

p(A 12) + I 2) • - ^p'(A) 

-Q 


and 


Q 

D ^ Jp(AIQ)-dA 

-Q 


(0.14) 


where D - limi^+oo {NIL) and K{A) = 0'{A)l{2n), viz- 


sinh(20 



2n sin(ju -i- i^) sin(/r - it^) 

for 

A > 1 


K(A) = . 

1 

7r(l -F 42) 

for 

A - 1 

(0.15) 


sin( 24 ') 

for 

-1 < A < 1 



In sinhiju + i^) sinh(// - i^) 


Yang and Yang did prove that the system of equations on the pair of unknowns (Q,p(A | Q)) does indeed 
admit a unique solution. They however did not prove the statement relative to the existence of the limit in the Ihs 
of (10.131) nor its value given by the rhs of (10.131) . 


Numerous works, starting from the pioneering handlings of Hulten lfT3l . did rely on the assumption that the 
limit of sums as in (10.121) exists and takes the form (10.131) . be it when in the Ihs there appear ground state Bethe 
roots or those describing a certain class of excited states above the ground state. In particular, such properties were 
used in the exact (but not rigorous in the sense introduced by Baxter |[T1) calculations leading to characterising 
the ground state energy and spectrum of excitations of the infinite volume XXZ chain |l5l[6l[T3l|25l|26l, testing 
the conformal structure of the spectrum of the XXZ chain l|3l[T6l[T7l|24l[33l[34l[35l, the algebraic Bethe Ansatz 
based calculations of the matrix elements of the reduced density matrix and, more generally, ground state 
correlation functions in this model ifTSl l20l or of the large-volume behaviour of the matrix elements of local spin 
operators taken between two excited states close to the ground state so as to name a few. Despite its 

importance due to the mentioned multiple applications, the existence and form of the limit (10.131) was not proven 
in its full generality so far. The only work which did address this question was the one of Dorlas and Samsonov 
HOl. The two authors focused on the case of the ground state Bethe roots and proved that indeed (10.131) does hold 
for -1 < A < 0 and for A > Aq > 1 where Aq was some explicit number. For -1 < A < 0 they could build 
their proof on a generalisation of the convexity arguments that were invoked by Yang and Yang relatively to the 
existence and uniqueness of solutions to (10.71) . Their argument was however limited to this regime since convexity 
does not hold anymore for A > 0. The two authors also managed to prove the statement for A > Aq > 1 by using 
the fixed point theorem for an auxiliary operator which was contracting for this range of the anisotropy. 

One of the results of the present paper is the proof of (10.131) for all values of A > -1 and for the class of 
real valued, particle-hole, solutions to the logarithmic Bethe equations (10.71) whose existence is established in 
Proposition 12.11 For simplicity, below, I only state the result in the case of the ground state Bethe roots. The 
general case can be found in Theorem 14. II 

Theorem. Let N,L ^ -i-oo in such a way that NIL —> D with D e]0; 1 /2]. Let q be the unique solution to (10.141) 
subordinate to D and {Tq}^ correspond to the set of Bethe roots parametrising the ground state of Then, 
given any bounded Lipschitz function f on R, it holds 

L 

^ ff(A)p(Alq)-dA. 

a=l 


(0.16) 






















In fact, this proposition is a corollary of a much stronger result established in the core of the paper : the 
existence of an all-order asymptotic expansion for the counting functions associated with such Bethe roots, see 
Theorems 13.51 and 13.41 The counting function contains all the fine details of the distribution of the da’s, so that 
obtaining this asymptotic expansion goes much further that the simple limiting result (I0.13I) . The idea of introdu¬ 
cing the counting function as a way to study the large-L behaviour of a given solution to the Bethe equations goes 
back to the work of l|9l|22l. The counting function formalism was further developed in the works iTTl l8l l2^ \2M 
what allowed to derive the first few terms of the large-L asymptotic expansion of the counting function associated 
with various configurations of Bethe roots associated with the XXZ chain. However, these derivations did build 
on various ad hoc hypothesis which, technically speaking, boil down to a densification property of the type (10.131) 
or a close variant thereof. In this paper I circumvent the use of such ad hoc hypothesis, hence bringing these for¬ 
mal asymptotic expansions to a rigorous level. It is important to stress that the method of proof introduced in the 
present paper neither relies on convexity arguments nor on fixed poinfs fheorem. The properfies I use are rafher 
general whaf makes fhe mefhod applicable, in principle, fo many ofher insfances of quanfum infegrable models. 

The paper is organised as follows. In Secfion [H I review some properties of solufions fo linear infegral equa¬ 
tions fhaf are relevanf fo fhe problem and esfablish certain auxiliary resulfs fhaf are of inferesf fo fhe analysis. In 
Secfion |2l I esfablish fhe solvabilify of fhe class of logarifhmic Befhe equafions of inferesf fo fhe problem. Then, 
in Secfion |3l I esfablish fhe main resulf of fhe paper, namely fhe exisfence of fhe asymptotic expansion of fhe 
counfing function. Secfion 0] is devofed fo fhe applications of fhe resulfs fo various problems fhaf arose earlier in 
fhe liferafure. In parficular, fhe densificafion proposifion is esfablished fhere. 

Appendix 13 confains an auxiliary resulf of inferesf fo fhe analysis. 


1 Some properties of solutions to linear integral equations 

Lef / c R be Lebesgue measurable. Lef Kj denofe fhe infegral operator on L^{J) acting as 

K7[/](d) = ^K(A-p)f(p)dp. (1.1) 

Here X(/l) corresponds to one of fhe fhree infegral kernels given in (10.151) . depending on fhe value of A. When 
A > 1, I will always assume fhaf diam(7) < tt, where diam(7) is fhe diamefer of J. Throughouf fhe paper, fhe 
dependence on A of fhe operators and infegral kernels will be kepf implicif since Ibis would nol bring more clarify 
fo discussion while weighfing down fhe infermediafe formulae. 

The purpose of fhis secfion is fo recall some known fads abouf fhe operafor id-tKy and fhe solufions fo specific 
infegral equafion driven by if. In parficular, I will discuss ifs invertibilify for any J and review several properties of 
solutions of linear infegral equafion driven by fhis operafor. Finally, I will prove an auxiliary resulf relative to fhe 
unique solvabilify of a non-linear problem driven by fhe operator Kj. This unique solvability will appear crucial in 
a later sfage of fhe analysis. 

Throughouf fhe paper, given a > 0, la sfands for fhe segmenf cenfred around 0 

Ia = V-a-,a]. ( 1 . 2 ) 


There are certain choices of fhe inferval J which make fhe operafor K j special, in fhaf fhe linear infegral equafions 
driven by id -i- K/ can be solved in a closed form by using Fourier fransforms or Fourier series. These choices 
correspond fo 7 = R when -1 < A < 1 and J = [-7r/2 ; 7r/2] when A > 1. If is convenienf fo infroduce a parameter 
L labelling fhe range of infegrafion /, corresponding fo fhose special cases : 


L 


oo if -1<A<1 
7r/2 if A > 1 


so fhaf 


4 - 


/oo - R if -1 < A < 1 

4/2 ^ [-:^/2 ; 7r/2] if A > 1 


(1.3) 
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It is a well known result that will be recalled in Lemma 11.31 that the operator id + Kj is invertible for all J. 
The resolvent operator associated with K/ will be denoted by Rj. It is defined as the operator on L^{J) such that 
(id - R/)(id + Kj) = id. The integral kernel of the resolvent will be denoted as Rj(A,ij). It solves the linear integral 
equation 

- y)Rj{y,^i)dv - K{A-fi) . (1.4) 


1.1 Some explicit solutions and their properties 

1.1.1 The density of Bethe roots 

The so-called density of Bethe roots is defined as the solution to the linear integral equations 
(id + K/J[p(* I q)] = ^p'(A) . (1.5) 

The density can be computed in closed form when the support of integration is 
Lemma 1.1. The solution to the linear integral equation 


takes the form 

Poo ( 4 ) — ^ 


1 


2 cosh (jtA) 


4p'(i) 

for 

-1 < A < 1 

Trw 

for 

A > 1 

for 

-1 < A < 1 

and 

A = 1 


(1.6) 


P7r/2{A) ^ 


neZ cosh [-(UTT - 4)] 


(1.7) 


The form of the solution, when -1 < A < 1 is readily obtained by taking the Fourier transform of the linear 
integral equation. When A > 1, one solves the linear integral equation by means of Fourier series expansions. This 
yields 


1 ^ 

~ 2^ 2 cosh(n^) ■ 

neZ ^ 

The expression (11.71) is obtained by applying the Poisson summation formula. These results appeared, for the first 
time, in ll^[J7l . It is clear from (11.71) that both poo and p;r /2 are strictly positive functions. 

1.1.2 The resolvent kernel 

It is readily seen that /?/,(4,p) only depends on the difference of its arguments. This integral kernel will be 
denoted below as R{A - p). 

Lemma 1.2. The solution R{X) to the equation (id + = K is an even function for any A > -1. 
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• For -1 < A < 1, /? admits the Fourier transform representation 


R(A) = 


-f 


sinh[(7r/2 - Ok]e 


-ikA 


dk 


cosh(^A:/2) sinh[(7r/2 - ^/2)k] An 


A = cos(^) with ^ £]0; n[ , 


(1.9) 


has the large A estimates R(A) = O^e and is a strictly positive functions when 0 < A < 1. 

For A = 1, /? admits the Fourier transform representation 


R{A) = 


-f 


e 2 . e 


—ikA 


dk 


cosh(^/2) 4jt 


( 1 . 10 ) 


has the large-A estimates R{A) = 0(2 ^) and is strictly positive on R. 

• For A > 1, /? has the Fourier series expansion 

I qIiiiA . g-|n|f 

RiA) = —\ - , A = cosh(0 with ^ > 0 , 

2n ^ coshfnn 

neZ 

is jT-periodic and strictly positive on R. 

The only non-trivial statement is the one relative to the signs of R. When A > 1, using (11.81) and 


( 1 . 11 ) 


PW 

2n 


z 

neZ 




nil 


7T 


one gets RiA) = | p„/ 2 iA - p)p'(ju)dp 

-nil 


( 1 . 12 ) 


from where striet positivity is manifest in virtue of (11.71) . When A = 1, one can recast the Fourier transform 
representation of R in the form R{A) = ^ ^ Poc,{p)v'{A-p)dp. The latter representation ensures strict positivity and 
readily yields the 0(2“^) estimates for the decay of R at infinity. Finally, after some calculations, when 0 < A < 1 
one recasts R as the convolution 


RiA) = 


n 


2 ^ 7^-0 


f 


sm 


TT-^ 


cosh ln(A - sinh pAtlFF}.] I’PizlPF] 


dp 


(1.13) 


TT-^ / \ TT-^ I 


which produces a manifestly strictly positive function. ■ 

The representation (11.131) was first found in llT/l while (11.121) is a straightforward application of the idea that 
leads to (Il‘l3l) . 

For the sake of further handlings, it is useful to introduce the integral operator L/ on L^iJ) defined as 


L/[/](2) - ^RiA-p)fip)dp. 


(1.14) 


If will be esfablished in Lemma [Ql fhaf id - L/ is inverfible. The resolvenf of fhis operator will be denoted £j and 
is defined by (id - L/)~^ = id + £y. 
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1.2 General considerations 

Lemma 1.3. Let 7 c R such that 0 < |7| < +oo if < A < 1 and diam(7) < n if A > 1. Then, the operators 
id + iC/ and id - Ly are both invertible. The integral kernel Rj^A, p) of the resolvent operator Rj satisfies to the 
bounds 

RiA-p) < Rj{A,p)<K(A-p)<0 for -1<A<0 (1.15) 

and 

K{A-p)> Rj{A,p)> R{A-p) >0 for A>0. (1.16) 

Furthermore, one has 2jc(A,p) = Rj(A,p), where = f\ J. 


Proof — 

For any g € L'^{J), one has the bounds 


ll^7[g]|lL»(R) < llgllL»(/) ■ sup r \K{p)\dp < 

AeK J 


J+A 


llgllL»(7) ■ sup ||i<:|lil(/) < 

/eJa 






where liJi /2 ^ [-|/|/2 ; |7|/2] and 

Ja = (/ c R : |/| - |7|) if - 1 < A < 1 and Ja - {/ c [-n-/2 ;n-/2] : |/| = |7|) when A > 1 . (1.18) 

The second bound is trivial for -1 < A < 1 and it holds for A > 1 since by 7r-periodicity of K one can reduce the 
integration along an interval of diameter at most n into one over an appropriate subinterval of [-7r/2 ;7r/2]. The 
bound holds since lA'I is even and decreasing on R"*”, resp. [0; 7r/2[, when -1 < A < 1, resp. A > 1. 

Now due to 


J K(ju)dp = 


if - 1 < A < 1 


and 


nil 

J K(p)dp = 1 

-nil 


if A > 1 


(1.19) 


the hypothesis on J ensure that II^IIl>(/|ji/ 2 ) ^ ^ Neumann series 


X n-i 

K{A-Vi)WK{Va 

n^L ’ a=l 


Va+l) ■ K(vn-p) ■ d"y 


( 1 . 20 ) 


for the resolvent kernel converges uniformly on R^. This establishes the invertibility of (id + K/). The one of 
(id - L/) is proven analogously using the properties of the resolvent R(A) established in Lemma fL^l 

When -1 < A < 0 the Neumann series for R/(A,p) is a sum of strictly negative terms. Hence the upper bounds 
given in (11.151) . For A > 0, following Yang and Yang lIT/l it is enough to observe that the integral equation for 
Rj(A, p) can be recast as 


(id + K4)[RX*,//)](2) - Kyc[7?X*,)U)](d) - K{A-p) i.e. (id - Lyc)[7?X*,)U)](d) = R{A-p) (1.21) 

where 1 remind that - f\ J. Since R{A) > 0 for A > 0, the Neumann series@ for the resolvent kernel £jc (T, p) 
consists of strictly positive terms what entails Rj{A,p) > R{A - p). The upper bound follows from K{A - p) > 0, 
Rj(A,p) > 0 and Rj(A,p) = K{A - p) - K 7 [R 7 (*,/r)](T). Finally, the equality 2jc{A,p) = Rj(A,p) between the 
integral kernels follows from the fact that A i-> 2jc(A,p) is the unique solution to the linear integral equation 
appearing to the right of (11.211) . ■ 


f. Its convergence follows from similar bounds and properties as those used for the one associated with K(A) 








Lemma 1.4. Let A > 0 and be two bounded and disjoint subsets of Let [J? 7 (+)] 7 (-) be the integral operator 
on acting as 


[i?7w]/-)[/]W - f RMA,p)f(p)-dp. 

j(-> 

Then, the operator ^id — [i?/+)]/-)) is invertible. Its resolvent operator [^/+)]/-) defined by 

^id - [i?y(+)]j(-)^^id + [.^7(+)]/-)) - id 
has a strictly positive resolvent kernel: \I%j{+)\j(-){A,p) > Ofor any A,p e R. 

Proof — 

For any g e it holds 


( 1 . 22 ) 


(1.23) 


||[R7w]/-)[^]|Lo„(r) ^ llgllL“(y(-)) • sup j |R/+)(d,/r)|d// 

/IeR kJ 

/-> 

^ llgllL“(/-)) • sup r \K{p)\dp < ||g||i<»(y(-)) • \\K\\[^i(i , ^ ) (1-24) 

ieR J y ) L tV(+)|/2l 

/+)+i 

where the second bound follows from the inequality (11.161) . The rest of the proof is carried out analogously to the 
one of Lemma [T31 ■ 

1.3 The magnetic Fermi boundaries 

When studying the observables of the XXZ chain in the thermodynamic limit NjL —> D, one naturally ends 
up with operators id + K/g where the endpoint Q is one of the unknowns of the problem. The endpoint q{D) which 
will be pertinent for describing the thermodynamics at given D and the so-called dressed momentum p{A \ q{D)) 
describing the momentum carried by an individual excitation over the model’s ground state arise as the solution 
to the below problem for two unknowns {Q, p{A \ Q)) : 

(id + K,g)[p(* I Q)](A) - ^ - £[0(4 - G) + 0(4 + Q)] and p(G I 0 = | . (1.25) 

It was shown by Yang and Yang llTTl that for any D e [0; 1 /2] there exists a unique q{D) e R'*' such that 
{q{D),p{A I q{D))) solves the problem given above. Furthermore, for A > 1 one even has q{D) £ [0;7r/2]. The 

cQ 

proof is relatively straightforward and builds on the fact that Q i-> J_g p'{A \ Q)dA is an increasing function of 
Q, c.f |[Tn for some details. The map D i-> q{D) is smooth and strictly increasing diffeomorphism from [0; 1/2] 
onto [0; -i-oo]. 

In the following, I will simply denote by q the endpoint of integration solving (11.251) . This endpoint will be 
referred to as the magnetic Fermi boundary. 

It turns out that in the intermediate steps of the analysis, it will become necessary to consider a slightly more 
general problem to (11.251) . namely one for three unknowns {Ql, QR,f{A \ Qi, Qr)) 

(id + ■Q,])[f{* I Ql, Gi?)](4) = ^ - £ • [0(4 - Qr) + 0{A - Qfi)] (1.26) 
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and 


KQr I Ql,Qr) = I KQl I Ql,Qr) - (1.27) 

with Qi^ < Qr, and the additional constraint \Qr - Qi\ < tt if A > 1. 

Proposition 1.1. Given any D G [0 ; 1 /2], the problem (ll.26l) - (ll.27l) admits the unique solution {q, -q, p(A \ q)), 
where {q, p{A \ q)) corresponds to the unique solution to the problem (11.25 1) . 

Proof — 

It is evident that {q, -q, p(A \ q)) provides one with a solution to the given problem. Hence, it remains to prove 
uniqueness. Let {Qi, QR,f{A \ Qr, Qr)) be a solution to (ll.26l) - (ll.27l) . It is convenient to distinguish between the 
case where both Qr and Qr are infinite, one of them is or both are bounded. Note that the first two situations can 
only arise when -I < A < I. 


Both endpoints are infinite 

First assume that Qr - -oo and Qr = +oo. Then, 

A 

f(A I -oo,+CX)) - p(A I +O0) = J' Poc(A)dA => p{±00 I +oo) = ±- 

0 

so that D = 1/1 and either this cannot be or one simply recovers the solution to the problem (11.251) . 


(1.28) 


One of the endpoints is infinite 

By symmetry, it is enough to deal with the case where Qr is bounded while Qr = +oo. For simplicity, denote 
simply the solution by f{A). An integration by parts ensures that (id + K[g^ ;+oo[)[/'] - pV27r. Then, it follows from 
the bounds (11.151) and (11.161) that 


< f'(A) < poc(A) for -1<A<0 and >/'(d) > poo(d) for 0<A<1. (1.29) 

Itt An 

These bounds ensure that / is bounded on R and strictly increasing. As such it admits a limit at ±oo. Hence 

|^(d-p)/(p)| < Ce^-^ for A<Qr so that ^hm^jK[e,;+oo][/](d)) = 0 (1.30) 

by dominated convergence. Thence 


lim /(d) = - 

/t—>—OO 2 . 


D TT-ni 




(131) 


Furthermore, due to the upper bound in (11.291 ) and lim^^+oo /(d) = D/2 it holds, for A> Qr 


/(d) - D/2 < Ce 


min ( 2 , j)a 


< Ce 


-A 


(1.32) 


so that 


K[GL;+oo][/-D/2](d) 


< C 


,-(A+/j) 


dp 


/i—»+oo 


0 . 


(1.33) 
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Thus recasting the integral equation for / in the form 


+00 

f(A) + fm- ^^)[m - I ] • = ^ - - 2l) ■ (1-34) 

Ql 

Taking the A +oo limit of the above equation yields and equation for D which implies that D = lf2. In its 
turn, due to (11.311) . this implies that lim^^_oo [/(d)] = -D/2 hence contradicting fiQi) - -D/2 since / is strictly 
increasing on 1. by (I1.29I) . 


Both endpoints are bounded 


In this case, it is convenient to introduce 

Qr - Ql Qr + Ql ~ 

Q ^ ^R ^^L ^ ^ ^ ^R ^^L ^ f{A + a\Qt,QR) ■ (1-35) 

One gets that f(A) satisfies 

(id + K/J[7] - - £ • [0(d -Q) + 0(A + Q)] and f(Q) - -f{-Q) - | . (1.36) 

The function / can be uniquely decomposed as /(d) = fi{A) + fp(A) where fp is even and /• is odd. It is readily 
seen that the functions fp and /• satisfy the linear integral equations 


(id 



p(d + a) + p(d - a) - D • [6/(d - Q) + 6»(d + Q)] 
p(d + a) - p(d - a) 


and are subject to the constrains 

fp{Q) - 0 and MQ) = | . 


(1.37) 


(1.38) 


The integral operator appearing above acts component-wise on the entries of the vector. Consider the equation 
satisfied by fhe even parf fp. If will be shown fhaf fhe consfrainf fp{Q) - 0 can only be satisfied if a = 0. Once fhis 
is esfablished, fhen (11.371) reduces fo (11.251) . whaf ensures uniqueness. In fhe course of doing so, one should freaf 
fhe fwo regimes -1 < A < 0 and A > 0 separately due fo fhe change in fhe sign of fhe integral kernel K{A - /i). 


-1 < A < 0 


By (11.151) . fhe infegral kernel of fhe resolvenf operator R/^ fo id -i- K/g is such fhaf /?/g(d,/r) < 0. Thus, due fo 
(11.371) . fp can be represenfed as 


/p(d) = 


p(d -I- a) - p(d - a) 
An 


- J /?/g(d,//)[p(/i 

-e 


da 

+ a) - p(ju - a) ■ — . 

An 


(1.39) 


Since 


A+a 

p(d + a)-p(A-a) = J p'Qj.) ■ d/r 
A—a 


(1.40) 
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has the same sign as a, it follows that, for any A £ R"*", 


if a>0 then fp{A) > —[p(A + a)-p{A-a)] > 0 

r • 

if a<0 then fp(A) < — [p(/l + a) - p(/l - a)] < 0 

Either of the two are inconsistent with the constraint fp{Q) - 0. Therefore, necessarily, a = 0 so that the problem 
reduces to (11.251) and on that account is uniquely solvable. 

• A>0 

We start by observing that the solution 0;;^, c.f (11.31) for the definition of l, to the linear integral equation 

A-¥CI 

(id + —~ takes the form (pt-aiA) = Jp,(ju) ■ dju . (1.42) 

A—a 

This can be seen as follows. When i - oo, the decay at infinity of the resolvent kernel R{A-p) and the boundedness 
of the driving term ensure that (poo-,a is bounded on R. When l = njl, it is readily inspected that <pni 2 -,a is 7r-periodic. 
Thus, for any l either because the boundary terms vanish (i = oo) or cancel out (i = njl), taking the derivative of 
the linear integral equation in (11.421) and then integrating by parts one gets a linear integral equation satisfied by 
0'.^ whose solution reads (p'^.J^A) = pi(A + a) - pi(A - a). The infegral represenfafion for (p^-^a then follows from fhe 
facl fhaf if is an even function. 

fp can be readily continued by means of fhe linear infegral equafion fo fhe real axis. Then, building on the 
identity 

(id + Yii^)[fp]{A) = (id + K/,)[^](T) - K/^f/pKT) with I^Q = h\lQ (1.43) 

it is easily seen, in virtue of (11.421) . that fp satisfies the linear integral equation 

(id - Lp^)[fp]iA) = ^(PM (1.44) 

where L/^ is as defined in (11.141) . As a consequence, if holds, 

IpW = \‘PM + f Ql^^(A,p)<Pr.a(M) • Y . (1.45) 

As a consequence, for any T £ R"*", 

if a > 0 fhen fp(A) > -<pi a(d) >0 for T > 0 

_ 1 • 
if a < 0 fhen fp{A) < -fi-aid) <0 for T > 0 

Again, either of the two are inconsistent with the constraint fp{Q) = 0 so that, necessarily, a = 0. ■ 
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1.4 Auxiliary functions 

In this last subsection I introduce a few auxiliary special functions that will arise in a later stage of the analysis. 
These functions are defined as solutions to linear integral equations driven by id + K/^ and thus depend on a free 
parameter Q. 

The dressed phase (f{A, | 2) is defined as the solution to the linear integral equation 

(id + I Q)]iA) = ■ (1-47) 

Above, * denotes the running argument of the function on which the integral operator acts. The dressed charge 
2(4 I Q) is defined as the solution to the linear integral equation 

(id + K/)[Z(H G)](4) = 1 . (1.48) 

The dressed energy £{A \ Q) is defined as the solution to the linear integral equation 

(id + K/)[£(* I 2)1(4) = e(4) with e(4) - /i - IJxk • p'(4) (1.49) 

where the constant xi^ depends on the anisotropy as 


Xa = sin^ 

for 

A = cos(0 

0 < ^ < n- 


Ta = 1 

for 

A = 1 


(1.50) 

II 

< 

for 

A = cosh(^) 

^>0 



Finally, the thermodynamic counting function is defined in terms of the dressed momentum p{A \ q) as 
^o(4k) = p(4|^) + | . (1.51) 

Lemma 1.5. The dressed phase is related to the dressed charge as 

^(4,212) - ^(4,-212) + 1 - Z(4k) and 1 + ^(2 I 2) - ^(-2 I 2) = (1-52) 

and also satisfies 

dAT{A,p I 2) - /?/e(4,/i) + Riq{A, Q)ip{Q,p \ Q) - Ri^(A,-Q)<p(-Q,p | 2) ■ (1-53) 

Also, the thermodynamic counting function satisfies ^q(4 | q) - p{A | ^) > 0 and thus is a strictly increasing 
diffeormorphism from 

R onto ^o(R \q) - [ - —(^ - D ); ^(^ - D) + d] (1.54) 

L TT 2 7T 2 J 

for -1 < A < 1 and from R onto R when A > 1. 

The proof of most statements is rather straightforward. In fact, the only non-trivial identity corresponds to the 
relationship with Z“^(2 I 2)- The latter was established by Korepin and Slavnov GTl and I refer to their paper for 
more details. 

For magnetic fields below fhe critical field 

{ 87cos^(^/2) -1 < A = cos(^) < 1 ^£[0;7r[ 

he = \ 2 (1-55) 

( 87cosh^(^/2) A = cosh(^) >1 ^ > 0 

the XXZ chain is an antiferromagnet. Furthermore, when A > 1, 0 < /i < h^P the model has a mass gap with 
h^^^ = 87sinh^(4'/2) whereas it is massless for h^^^ < h < he. 
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Proposition 1.2. Assume that the magnetic field h satisfies to the bounds 

0 < h < h(. for - 1 < A < 1 and h^^^ < h < he for A > 1 . (1.56) 

Then, there exists a unique solution Qp, called the Fermi boundary, such that e{Qf \ Qp) = 0. Furthermore, it 
holds 


siQ I G) < 0 for Qe]-Qp- Qp[ and s{Q | G) > 0 forQel] 


Qf ■ 


(1.57) 


This result has been established in ifTTI when -1 < A < 1 and the technique of proof can be readily extended 
to the regime A > 1.1 refer the interested reader to that paper for more details. 


2 Solvability and boundedness 

2.1 Existence of solutions to the logarithmic Bethe equations 

Proposition 2.1. 

• -1 < A < 1 

Let ^ G [0; 7r[ parametrise A = cos(^) and 


hi <■■■ <h„ , hael\-,NJ 


Pi <■■■ < Pn , Pa ez\|[l; , 


be ordered integers such that 

- Hfll) 

n L L 

Define the set of integers {(a]i by 


Pi - 1 

L n ^2 


N-\ 

L 


) 


and 


n ^2 


Ns n 

t)^j- 


ta=a for ae'i\\ N'^\{hi,...,hn] and h,, = Pa for a=\,...n . 


(2.1) 


( 2 . 2 ) 


(2.3) 


Then, the system of logarithmic Bethe equations (10.71) with p and 6 as defined in (10.81) or (1091) admits a solution 
such that all Bethe roots {Ta}^ are real. Furthermore, for ^ G [7r/2 ; 7r[, this solution is unique. 

• A > 1 

Let ^ g] 0 ; +oo[ parametrise A = cosh(^) and ha e E 1; and pa € be any ordered integers, viz. 

hi < ■ ■ ■ <hn and pi < ■ • ■ < p„. Then, the system of logarithmic Bethe equations (10.71) with p and 6 as defined in 
(10.101) admits a solution such that all Bethe roots {Ta}^ are real. 

Note that the proposition only stipulates the existence of solutions in the case of general A > -1. Uniqueness 
only holds, a priori, for 0 > A > -1. Furthermore, when A > 1, the proposition only states that the da’s are 
real and says nothing about the domain to which they belong. The main peculiarity of the A > 1 regime is that 
two solutions of the logarithmic Bethe equations which differ by translations of n, namely Aa = A'a + man for 
a = I,... ,N and some G Z define equivalent solufion fo fhe Befhe equations. Hence, disfincf sefs of infegers £a 
and {'a do nof necessarily allow one, when A > 1, fo disfinguish befween inequivalenf solufions. So as fo deal wifh 
only one represenfafive one should, in fine only focus on fhe solufions belonging fo some fixed inferval of lengfh 
n, say ] - 7r/2 ; 7r/2]. This is, however, somefhing fhaf should be done a posteriori, affer having builf fhe solutions. 

The proof of fhis sfafemenf basically follows Yang-Yang’s argumenf for fhe ground sfafe Befhe roofs when 
-1 < A < 0. The idea consisf in infroducing a funefion p i-> S lip) on such fhaf fhe logarifhmic Befhe Ansafz 
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equations correspond to the conditions for a local minimum of S t- The main new observation introduced here is 
that even though the function S (is not convex for A > 0, it still blows up at oo provided that (12.21) holds, and hence 
admits a minimum. 

Proof — 

• -1 < A < 1 

Following the reasoning of Yang and Yang, one introduces the function 


N or \ . N 

a=l a,b=l 


defined in terms of 


1 w ^ n 

— ^ &Qia-Pb) - with na = 4 

a.b=l a=\ 


A^+ 1 


(2.4) 


A 

-f 


P(A)= p(p)dp = (^-4IT| + - 


0 ( 1 ) 

0(ln|d|) 


A 


and 0(d) = 0(ju)dju = (7 t-20\A\ + - 


0 ( 1 ) 

0(ln|d|) 


■ (2.5) 


The remainders appearing above are such that the first line holds pointwise in ^ €]0; 7r[ while the second line 
holds for ^ = 0. This convention will be carried on until the end of this proof. 

It is easy to see that the logarithmic Bethe equations (10.71) appear as conditions for the existence of a critical 
point d of 5 ((ji) : 


4- • Sifx) 

\i=ii 


N 


P(dfl) 1 V t 

h=\ 


( 2 . 6 ) 


Thus, it is enough to show that Sdji) —> +oo as > +oo so as to ensure the existence of a minimum of S dji) and 
hence the existence of a solution to the logarithmic Bethe equations. 

It appears convenient to study the behaviour of S e(ji) when fi goes to infinity along a ray. In such a situation, 
there exists cr € such that = (fio-n), ■ ■ ■,Pct(N)) goes to infinity as 


Ha-^tv with V = ( vi,..'^.,vi , ..., y„... ,vj vi < • ■ • < V;-< 0 < V;.+i < ■ • ■ < v, . (2.7) 

fTi terms as terms 


There t —> +oo and v e S^, the N-dimensional sphere. Given n as above, when t —> +oo, it holds 

n - ( ^ 7r-2(('^ - j 0(1) 

+ Zj o!(V(-t-— 2^aja({v(-vj) - -)^V(ne + 


S({p) = 


t=\ 


£=r+\ 


j<t 


e=\ 


O(ln0 


■ ( 2 . 8 ) 


Here, we have set 


at t 

nt= ^ ncr(;) and ar = ^ , oq ^ 0 . 

j=at-i + l p=\ 

Therefore, using that 

5 S 5—1 

^ ajadve -Vj) = ^ atvtat-\ - ^ ajVj (77 -af) , 

j<t t=i j=\ 


(2.9) 


( 2 . 10 ) 
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one can recast the large t asymptotics of S lijx) as 


= t {(vr - vm ) • (+ v, • ( ^ 
e=\ ;=i y=i 

S S 

+ Vf_i) • (Y ^'■+1 

i=2+r j=£ 

where 

^-4' A^(n- - 20 _ a^_i +af 7r-2^ | _ ^ 

27r 27rL 2L n ! L 



(z -r‘)) 


j=l+r 


0 ( 1 ) 

0(ln t) 


One will have S i(ji) + ex? for any ^ co provided that 


( 2 . 11 ) 


( 2 . 12 ) 


{ S 

Yj Y'’ ^ ® for all ^ e H 1; rj and ^ >0 for all ^ £ |[r + 1; 5 ]] (2.13) 

;=i j=e 


this for all choices of point v e S^, viz- for all admissible values of r, s and of the a^’s. These equation impose 
constraints on the integers ria- Indeed, using that 


{ 

p=\ 

one gets 


S S 

YjCtpif^p-i + <^p) - IN'ae-'a^ with He = Yj<^k = N-a(-i (2.14) 

p={ j=e 


V (-1 = { ^ + A^(7r - 2^ x 2 (^ ~ 2^ _ ^ kio-(k) 

j 27r 2n'L ' ^ 27rL Y L 

]=\ k=\ 

and also 


(2.15) 


V (+) ~ ~ ^ N{n-2^)^ _2-^(^“2^) ^ na-{N-k+\) 

j=i k=\. 

Now by running through all the possible types of inequalities (12.131) . one concludes that S lipi) 
^ > 00 , provided that, for any c H 1; NJ, the integers satisfy 


(2.16) 


+00 for any 


m 


( 


n-i^ N{n-10\ 2^-2^ 
In 2nL ' 2 tiL 


< 



< m 


( 


2n 


N{n-20 . 

2nL > 


+ m 


Tk-2i^ 

2nL 


(2.17) 


with m = #J'. 

Since, when (12.171) are satisfied, S((p) +00 for any // —> 00 , it follows that Se admits at least one minimum 
on at some point A. The coordinates of the point realising this minimum satisfy to the logarithmic Bethe 
equations, hence ensuring the existence of solutions. 

It now remains to check that the constraints (12.171) are always satisfied provided that the bounds (12.21) hold. 
Let riQ be such that 


Pi < ■ ■ ■ < Pno <0 < N < p„g+i <■■■ <pn - 


(2.18) 


16 

























Then, for any c [[ 1; A/” J, #J' = m, agreeing upon rriQ = min(/n, no), it holds 


mo 


> Y^{Pa- 

aej 


A^+ 1 


m—mo 


a=\ 


a=\ 


+ ^ -'« + 'wo-1) + 


m{m + 1) m(A^ + 1) 


(2.19) 


where the lowest bound has been obtained by using that Pa > P\ + a - \.\t thus follows that the lowest bound in 
(12.171) holds provided that 


-m- 


n - ^ N m 


(I i\ mx niQ 

+ < —{pi-m + mQ-\) 


( 2 . 20 ) 


Since mo < nr for m G K 1; no 3> these bounds are clearly satisfied as soon as (12.21) holds. Further, for m > no, 
the bounds (12.201) will hold provided that they hold for the smallest possible choice of pi compatible with (12.201) . 
namely by taking it equal to the Ihs of (12.201) when at m = mo = 1. This translates itself on the condition for no : 


, jTT-^^l 77 m + no^^ no^ ^ .n 

(m + _^) - -J - - l)_ < 0 . 


L 


ttL 


( 2 . 21 ) 


Since m > no, for the above to hold, it is enough that the latter always holds for 0 < ^ < n/l and imposes the 
constraint 


TT-^/l N 


L> 


lL(i - ley'll 
- Tr\2 T > ~ L 


( 2 . 22 ) 


which is clearly satisfied provided fhaf n is bounded according fo (12.21) . One can repeaf fhe same reasoning relati¬ 
vely fo fhe upper bound in (12.171) whaf leads fo fhe sufhcienf consfrainfs 


-m- 


TT — {/I N m\ 


/ r rv m\ mt 

- -p -)- y) > - N + m- m-i-j wifh ni'^ = mm(m, n - no) . 


For fhe laffer fo hold, if is enough fo have 


TT-C/l N !■ 


/I IS i\ mt n-L /L IS \ n 

-riPn-N + m.-m) and ( )> 


/I _ 

2^-7r''2 L' 


n 

L 


(2.23) 


(2.24) 


Clearly, bofh bounds hold if (12.21) holds. 

Now, regarding fo uniqueness, observe fhaf for ^ G [7r/2; S e(ji) is sfricfly convex since if has a sfricfly 
positive defined Hessian mafrix 




Se(M) - -jY,K(pk-Pt)]Sjk + 


N 


dpkdp 


K(pj- Pk) 


In 


(2.25) 


} “ b=\ 

Indeed, given (hi,., his) ^ one has fhaf 

N q2 N 


i,k=\ 


j 


k=\ 


: V'(Pk) 

2n 


N 

Z 

. 74=1 


(hj-hk)^ K(pj-pk) 


> 0 . 


(2.26) 


Thus S t(p) admits a most a single minimum, in this range of ^’s. 
• A > 1 
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The analysis follows basically the same lines as for -1 < A < 1. Defining Ar as in (12.41) with, now, the 
functions P and 0 being given by 


A ^ A ^ 

P(A) - I p(M)dju = ^ + 0(|T|) and 0(d) = J = ^ + 0(|d|) 


A 


(2.27) 


one gets that the logarithmic Bethe equations do arise as necessary conditions for a local extremum of S (. Sending 
/I to oo exactly as it was done for -1 < A < 1 leads to the asymptotic behaviour 


■2 ^ 


^b) + + ^^bf] + 0(0 . 


a,b=\ 




(2.28) 


Since LjN >2, this ensures that S i blows up at infinity and thus admits at least one minimum. 


2.2 Boundedness of solutions to the logarithmic Bethe equations 

I now prove that for -1 < A < 1 at density D e [0; 1 /2[ and under certain restrictions on the allowed range for 
the Pa’s and n, the solutions to the logarithmic Bethe equations are bounded uniformly in N,L. Boundedness also 
holds for A > 0 and D e [0; 1/2] under slightly different restrictions on the pa’s. When D = \j2 and -1 < A < 1, 
the Bethe roots are not bounded any more. I establish bounds on the proportion of roots lying outside of a compact 
of large size. 


Proposition 2.2. 

• Let -1 < A < 1 

Let N,L be such that N/L D e [0;l/2[. Parametrise A = cos(/i') with ^ e [0;7r[ and let 
correspond to a solution of the logarithmic Bethe Ansatz equations (lO) where the integers p\ < ■ ■ ■ < p„ are 
subject to the constraint 


^ ~ h > Pn-lV Pi - 1 ^ ^ - ^ /1 _ 1 \ 

7r''2 L L ’ L ^ Tx ^2 L > 


(2.29) 


and n is fixed and N independent. Then, there exists a N,L independent constants A > 0 such that 


\Aa\ < A for any a £ [[ 1; A^]] such that f’a £ II 1; 3 uniformly in (N,L) . 


(2.30) 


• Let A > 1. 

Let N, L be such that NIL —> D £ [0; 1/2] and ei Ti,..., div correspond to the solution of the logarithmic 
Bethe Ansatz equations (lOTl) where the integers pi < • • • < p„ are such that n is fixed, L independent, and 
|pa/A| < C for some L-independent C > 0. 

Then, there exists a N,L independent constants A > 0 such that 

|da| < A for any a £ (1,..., A) uniformly in (A, L) . (2.31) 

Note that given the constraints on the Pa’s, solutions to the logarithmic Bethe Ansatz equations do exist in 
virtue of Proposition 12.11 

Proof — 
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The proof splits in three steps since one has to distinguish between -1<A<0, 0<A<1 and A > 1. The 
distinction between the first two regimes is due to the change in the sign of the kernel K{A). The last regime has 
to be treated separately due to the change in the periodicity properties of the involved functions. 

• A > 1 

Let be the reordered Tq’s, namelyS {da}^ = {da}^ and di < • • • < d^?. Also, let {ia] denote the corres¬ 
ponding reordering of the f’^’s. Let 


ad) = — p(d) 
2n 


1 ^ 


2nL 


Aa) + 


a=l 


N+l 

2L 


(2.32) 


be the counting function built up from this solution to the Bethe Ansatz equations. 
Since the d^’s are real, by 7r-periodicity of the functions, it follows that 


P(d)| < :^P'(0) + V(0) . 


2n " ' 2 

Let p be such that {p + l)7r > d/^ - di > pn. Then, since ^{x + kn) - ^(x) = k{L - N)/L, one has 


(2.33) 


In - 


= ad;v)-adi) - adi+p7r)-adi) + j r(d)-dd > p 

A\+pn 


■/?'< 


~ ^p'(O) - . (2.34) 


This yields the upper bound on p since the ^^’s are bounded in L. Thus there exists an L-independent constant 
C > 0 such that C > |dAr - di | > 0. Since 6 is increasing on K., it holds 


1 . 1 N+\ N ~ C 

—V{^n) < 7max{4}- — + TT-pOiC ) = — 

2n L 2L 2nL 2n 


In < (P)‘^(C') - A 


(2.35) 


The lower bound on A\ is obtained analogously. 

• A = cos(a with ^ £ [0; 7r/2] 

Reorganise the Bethe roots as 

{dfl}f ^ Klf with Va = Aq^ and v„+i < • ■ • < va? with 4^e[[l;A^]l for a ^ n-i-1, ■ • ■, ■ (2.36) 


Since 9 is increasing in this region of ^’s and is bounded by tt - 2^, one has 


P(vjv) 

2n 


< + 


A- 1 
2L 


2nL 


N . 


Since iq^ < N 


v{vn) n-p N n-ip 


2n 


< -^ • -r < - -{D + e)^AN< (p)-V2(;r - + f)) ■ 

n L n ' ' 


where e > 0 is such that A/L <D + e < 1/2 and L is taken large enough. 
The bound on di is obtained analogously, leading to 


P(v«+i) 

2t: 


> 




(D -I- e) o Vn+\ > -p '(2(7r - ^(D -i- e)) . 


f. There cannot be two To’s that are equal since this would contradict that the L’s are pairwise distinct. 


(2.37) 


(2.38) 


(2.39) 
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Thus, the claim holds with A = p“^(2(7r - ^){D + e)). 

• A = cos(^) with ^ £ [7r/2 ; 7r[ 

It is a classical fact that if {Tq}^ solve the system of logarithmic Bethe equations associated with the choice 
of integers £a, then the set of roots {Ha)^ with Ha - solve the system of logarithmic Bethe equations 

associated with the choice of integers Therefore, the claim will follow as soon as one shows that 

limsup^^^^[y/v] < +oo, where the Vq’s are as defined in (I2.36I) . 

The proof goes by contradiction. Hence, assume that limsup^^^^[y/v] = +oo. In other words, there exists 
a sequence (L„,,Am), Nm,Lm —> +“, Nm/Lm —> D and the associated Bethe roots associated with the 

subsequences and chains of L,„ sites such that lim„,_j+oo = +oo. Fix Q> 0 and set 


e E” +1; I : > e) ■ 

i V f'f'i fl 

Define further 

bo = lim sup b^™^ and b = lim sup bn . 

m-^+oo ^—>+00 


(2.40) 


(2.41) 


By construction, one has 1 > b > 0. Furthermore, by definition, there exists an increasing sequence Q{, 
Qi —> + 00 , such that bg^ —> b. 

i—^+OO 

Assume that 1 > b > 0 and pick any e > 0 such that b - e > 0. Let {q be such that 


y£>£o l'-|<bQ, <b + |. (2.42) 

Also, given a fixed £ > £q, there exists a sequence r i-> mr{£) and ro(f’) e N such that for any r > ro(f’) one has 

+ | (2-43) 


One can slightly improve the rhs of the bound. Indeed, since bg^^ < b + e/2, there exists mo such that for any 
m > mo, one has b^™^ < b + e. As mr{£) goes to infinity with r, there exists 

r[{£) such that Vr > r[{£) one has m^(f’) > mo . (2.44) 

Thence, setting ri(£) = max(ro(£), r'^i£)), one has for any r > ri{£) that 

(b + e){NmAf) -n) > #{a £ 11 + n ; I : j _ (2.45) 

Since, the total number of the T’s is the above inequality implies that 

#|a £ H 1 + n ; ]] : y^^^O) < . max|0,1 - e - b) . (2.46) 


a ^ [A^«,,(r)(l + £ - b)] + 1 and a = maxjO, - e - 6)]) 

where [.] stands for the integer part. Then, 

^£q,, < -{Nm,(i) -a+ l){Nmr({) - a - n + \) . 

k=Ti k=a+n 


(2.47) 


(2.48) 
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Therefore, after summing up the logarithmic Bethe equations involving vj,..., invoking the above bounds 

and using that 6 is odd, one gets 


{NmAi) - a + \) 


a — \ 


N„,r(e) 


NmriO a-l 


2L 


mAf) 


> — Y ^— y y e{ 

27r Zj V /c ) InLmAt) ^ ^ ^ 


1 


In 

k=a 
NmrV) 9. 


^(mAt)) _ yM))\ 


k=a p=a+\ 


2nL 




(£)) JmA£)) 


•mAt) t - „ 1 

' k=a p=i 


) 


Nmr{t) 


2nL, 


miAi) 


k=a p=l 


,{mrm (Mrii)) 


)■ (2.49) 


Using that is strictly increasing so that for ^ > p > n one has 

(12:4^ 


,{mAt)) _ AmAt)) 

/, y n 


> -0(0) = 0 and that by 


\> Qe fox ke la ■, NmAC) 3 and At < Qe^ for ^ e [[ 1; qj 
by (12.451) . one arrives to the bound 

a - 1 ^ ViQt) q ^ ^ 2^-n 

T7- - - ~ - T~r -” • 

ALmAt) An InLmAt) AnL^At) 

Taking the r —> +oo limit, it follows from 


(2.50) 

(2.51) 


lim {alLmM)) ^ D ■ (1 + e - b) and lim {qlLmAt)) ^ D • max {0, (1 - e - b)) , (2.52) 

r^+oo r—>+00 

that 

|(1 + e - b) > - £ max |0, (1 - 6 - b)) me - Qe,) • (2-53) 

Sending first { +oo and then e —> O'*' leads to 


D- 


1 


> 




D 


n-2C 


(l-b) 


VIZ. 


D(1 


b) > - . 
2 


(2.54) 


2 ~ 2n 2n 

However, the last inequality cannot hold since, by hypothesis, 0 < D < 1/2. 

Thus, one necessarily has 0 = 0. This however does not yet guarantee that the roots are bounded from above 


since a small portion thereof can escape to +cx3. Thus assume that one has A 


logarithmic Bethe equations for A 


(mAt)) 

Nmr(e) ’ 


(mAO) 

'N,nr(,e) 


and using similar bounds one gets 


1 ,N, 


mA£) 


u 


■^mA£) 


P(^iV„,Ko) 

2n 


2nh, 


•mAA 




n(2^ - n) 


2nL, 


•mAt) 


+ 00 . Then, staring from the 


(2.55) 


with a and other quantities as defined above wifh fhe difference fhaf, now, b = 0. One can send r —> +cx3 on fhe 
level of fhese bounds and fhen relax e ^ 0 leading fo (12.541) wifh b = 0. This yields fhe soughf confradicfion. ■ 
In fhe resf of fhis subsection, I focus on fhe case -1 < A < 1 when D = 1/2 and obfain bounds on fhe 
proportion of Befhe roofs for fhe ground sfafe which lie away from some segmenf 1\. More precisely, one has fhe 

Proposition 2.3. Let -1 < A < 1 and NIL = D = 1/2. Let {Tq}^ be a solution to the logarithmic Bethe Ansatz 
equations associated with the ground state choice of integers ta - with a £ |[ 1; Then, for any e > 0 there 
exists Af, Lq > 0, such that for any L > Lq, 


ca < e 


where ca = —#|a £[[1;A3 : /laeR \ /a} 


corresponds to the fraction of Bethe roots lying away from the segment I\. 


(2.56) 


Proof — 
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• 0 < A < 1 


Straightforward bounds analogous to those developed in the proof of Proposition I2.2l lead to the upper bounds 


N-a 1 
L ^ H 


+ 00 



and 



■oo 


Assume that Aa > A. Then, one gets the upper bound 


N-a 

- < 

L 


+0O 



A 


and setting k\ 


+ CO 




+ 1 


(2.57) 


(2.58) 


where [*] denotes the integer part, it is clear that there can be at most distinct integers £ H 1; ]] satisfying 

to this constraint. Likewise, if Aa < -A, one has 


a 

L 



, du 1 


Thus setting 


.jp'W-l 


+ 2 


(2.59) 


it becomes clear that there can be at most distinct integers la £ E 1; A ]] satisfying to the above constraint. 
Hence, all-in-all, one gets the bound 



R\/a 



(2.60) 


The conclusion then follows from the fact that the rhs of the inequality approaches 0 when L,A^ -i-oo. 


• -1 < A < 0 

If the sequence {Aa\^ is bounded uniformly in A, then the statement holds simply by taking A^ > M, where 
M is a bound on the magnitude of the Bethe roots. Else, one reasons as in the proof of Proposition 12.21 so as to 
establish that b = 0 with b as defined in (12.411) . One gets similar bounds relatively to the proportion of roots lower 
that some fixed paramefer Q. The resf is sfraighfforward. ■ 


3 Asymptotic expansion of the counting function 

3.1 Leading asymptotics of the counting function 


Given a solufion {Aa\^ to the logarithmic Bethe equations associated with the choice of integers it 


N 


IS 


convenient to introduce the associated counting function 


ad) - — p(d) 
2n 




< 3=1 


N + 1 
2L 


(3.1) 


By construction, this function is such that ^(Aa) = 7’^ for a = 1,..., A. The purpose of this section is to provide 
an alternative to (13.11) characterisation of the counting function valid in the large-L regime. This characterisation 
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is obtained by means of a non-linear integral equation. The very idea goes back to the works of De Vega, Woyna- 
rowich |'9] and Batchelor, Kliimper |[22]| and was further developed in the works 1711231. The non-linear integral 
equations obtained in the earlier literature were obtained by assuming that the counting function satisfies certain 
properties such as being strictly increasing on the real axis and, on compact subsets of R, having a bounded from 
below derivative, ^' > x > 0, this unformly in L. The main point of the method is that once a non-linear inte¬ 
gral equation is taken for granted just as certain amount of properties of its solution, then it is relatively easy to 
compute, order-by-order, the coefficients of its large-L asymptotic expansion. The assumptions which allow one 
to derive the non-linear integral equation for the counting function and which also allow one to derive its large-L 
asymptotic expansion could, in the best case scenario, be verified a posteriori, namely on fhe level of fhe obfained 
form for fhe large-L asympfofic expansion. This only allowed for a consistency fesf of fhe calculations. 

The main inpuf of fhe analysis fhaf I develop below is fo sef techniques allowing one 

i) fo prove fhaf, for L-large enough, fhe counting function can indeed be characterised as fhe unique solution 
fo a non-linear infegral and fhaf if does indeed satisfy fo fhe expected properties, in particular, fhaf if is 
sfricfly increasing on R.; 

ii) fo demonsfrafe fhaf fhe counfing funcfion admifs a large-L asympfolic expansion up fo o(L“^) corrections. 
In ofher words, fhe framework developed below allows one fo step ouf of fhe formal handling of asymptotic 
expansions. Prior fo sfafing fhe resulf and going info fhe defails of fhe proof, I need fo infroduce several building 
blocks of fhe non-linear integral equafions salisfied by 

Given D = V/L € [0; 1/2], here and in fhe following will denote fhe unique solufion fo fhe magnetic Fermi 
boundary problem (11.251) associafed wifh D. In ifs furn, q will denofe ifs fhermodynamic limif, viz- the unique 
solution to the Fermi boundary problem (11.251) associated with D - lim D. 

Definition 3.1. Let f be real analytic function such that it that is a biholomorphism on some open neighbourhood 
of a segment ofM. and such that its range contains, for some ao > 0 small enough, the closed set G{ao) depicted 
on Fig. [7] Then, one can construct the three below non-linear integral operators 



(3.2) 


and agreeing upon 


qAH =f , q^[f] =f 


(3.3) 


one has 


?«[/] 




and 


?l[/] 


(3.4) 


These operators, build up a "master” non-linear integral operator as 


3 



(3.5) 


a=\ 
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Figure 1 - Closed set G(a) delimited by its boundary r+ U r_. 


Just as in Section [H I will consider 

hi < ■ ■ ■ < h„ , ha e l[l; NJ and pi < ■ ■ ■ < , pa e Z\\[l ; N]\ 


(3.6) 


where n £ N is A/^-independent. However, the integers ha and pa can depend on N provided that they satisfy to 
certain bounds : 

• for-l<A<l,A = cosC^"), in addition, the integers {pa}” are assumed to satisfy to the additional constraint 


TT-^.l N-U Pn-N Pl-l TT-^/l N-\ 






n L L ’ L ^ 7T ^2 L 

• For A > 1, the integers (pa)” are solely assumed to be bounded as 

|pa/A| < C for some L - independent C > 0. 

In the course of the proof, it will be useful to introduce the sets 

S^{J) = (z £ C : |3(z)| < 77 , %{Z) £ /) 


and 


^ (zeC : |3(z)| < p, d{%{z),J) < e} 


(3.7) 


(3.8) 


(3.9) 


(3.10) 


where / c R is a subset of R and d is the distance between subsets of R induced by the Euclidian distance. 
Finally, it will also appear convenient to introduce the parameter defined by 

-1<A<1 = ^lA A = cos(^) 0 < ^ < n- 

A-1 = 1/4 . (3.11) 

A>1 = ^1A A = cosh(^) ^ > 0 
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Theorem 3.2. Given A > -1, let n eN be fixed and the integers {{a\^ be defined in terms of integers {paY\ 
{ha}1 according to (I2.3I) . Assume that these integers satisfy, depending on the value of A, the conditions (I3.7I) - (I3.8I) . 
Let N, L be such that 0 < D = NIL < 1/2 and go to infinity so that D D with 


D e [0; l/2[ for - 1 < A < 1 and D £ [0; 1/2] for A > \ . 


(3.12) 


Finally, let {/Iq}^ denote a solution to the logarithmic Bethe equations subordinate to the associated integers {€a]^- 
Then 

i) there exists Lq large enough such that, for any L > Lq the associated counting function (ITT]) is a strictly 
increasing diffeomorphism from 


R. onto 


TT-^/l N-. \ N 1 n-qi /7A 

TT ^2 ~T''^2L’T^2L^ jt ^2 ~ Tl 


(3.13) 


when -1 < A < 1 while, for A > ^ is a diffeomorphism from R. onto R. 

ii) There exists an open neighbourhood U of[-q ; q\ and an L-independent open neighbourhood Vd o/[0; D] 
containing G{a) delimited by the curves r+, c.f Fig. [7] for some oq > 0 such that ^ : U i—> Vd is a 
biholomorphism. 

Hi) There exists C > 0 such that 


ll^(*) -^o(* I - Y k) = k) + y ■ (3.14) 

Above, * denotes the running variable of the functions, 
iv) The counting function solves the non-linear integral equation 

JiA) = fo(4 I?) + ^ I {7pJl ; {7hJl) + 91^[?](4) , (3.15) 

with 

ae{hu...,h„,pi,-..,Pn} ■ (3.16) 

The function OA is expressed in terms of the dressed phase and charge as 


^ m n 

I [yfif ; [zaffi = 2^^ + I - Z I ^ + Z I ^ 

a=l fl=l 

and the operator 51a? is as given in Definition 13. il 

v) The non-linear integral equation (13.151) admits a unique solution in the class of functions satisfying to ii) 
and Hi) when A> 1 or-1<A<1 and, on top of the previous conditions, the Pa jL all belong to an 
L-independent compact subset of 


n 



\D + 


n 



(3.18) 


The constant C in (13.141) and Lq appearing above only depend on n and are uniform in D belonging to compact 
subsets o/[0; l/2[/or -1 < A < 1 and throughout the segment [0; 1/2] for A > 1. 
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In virtue of Proposition 12.1[ given integers {/Jq}” and {pa}” satisfying to the constraints (I3.7I) - (I3.8I) . depending 
on the value of A, there are always solutions to the logarithmic Bethe Ansatz equations. When -1 < A < 0, 
the solution was shown to be unique but, a priori, when A > 0, there could exist more than one solution to the 
logarithmic Bethe Ansatz equations associated with this given choice of integers. The statement of the theorem 
does hold for any such solution. 

Proof — 

For each value of N,L, one is provided with integers {haY{ and {pa}”, which, as mentioned, in most cases do 
depend on N and L. In virtue of Proposition 12. 1 1 these give rise to a sequence {Tq}^, indexed by L, of solutions 
to the logarithmic Bethe Ansatz equations. By Proposition 12.21 equation (12.311) . there exists and L-independent 
A > 0 such that : 

• Idfll < A for a £ H 1; if A > 1, 

• \/ia\ < A for a £ H 1; such that L^{Aa) £ E 1; in the regime -1 < A < 1, 

this uniformly in N, L. This sequence of Bethe roots {Aa}^ defines the associated counting function through (13.11) . 
therefore giving rise to a sequence (in respect to L) of counting functions f. These counting functions are readily 
seen to be holomorphic functions on the strip Furthermore, straightforward bounds show there exists an 

L-independent constant B > 0 such that 

(3.19) 

Since ^ is a sequence of holomorphic functions on Sxjyhh) that are uniformly bounded in L, by Montel’s theorem, 
it admits a converging subsequence which converges, in the sup norm topology on compacts subsets, to a 
holomorphic function on 

The strategy of the proof consists in characterising the limit of such a convergent subsequence. It will be 
shown that ^g necessarily coincides with the so-called thermodynamic counting function ^o(' I ^)^ c.f. (13.141) . Wor¬ 
ded differently, any converging subsequence of ^ has the same limit. Since, by Montel’s theorem, any subsequence 
of ^ admits a converging subsequence, f necessarily converges to \ q). Once that the convergence is establi¬ 
shed, the form of the non-linear integral equation (13.151) satisfied by ^ follows rather easily from the properties of 
the limit fo. A straightforward investigation of the non-linear integral equation ensures the uniqueness of solutions 
for L large-enough and the bounds (13.141) . 

In order to lighten the notations, I will subsequently drop the subscript e in all the considerations namely the 
subsequence and its limit will still be denoted by ^ and f. Since is holomorphic on hh, it admits, a finite number 
of zeroes and thus may only change signs a finite number of times. 

In the first part of the proof, I will assume that ^' > 0 and show that this property is enough so as to characterise 
the limit. In the second part of the proof, I will establish the uniqueness of solutions to the non-linear integral 
equation. Finally, in the third part of the proof, I will rule out the possibility that f has zeroes or changes sign on 
[-A;A]. 


• > 0 on [-A ; A] 

In virtue of Proposition lA.ll there exists q,e > 0 such that 
^ ^ fiS2,i,2e{lA)) is a biholomorphism Satisfying ^(52;,,26 (/a) n H*) 

Furthermore, for L large enough, it holds that ^iS 2 r^, 2 f;{lA)) ^ d ^(/a) and 

? : <S2,,2.(/A)n?-'(^(5,,g(/A))) ^ ^{SnA^A)) 


^(‘S2,,2.(/A))nH± . 


(3.20) 
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is a biholomorphism. Besides ^ is strictly increasing on /a+e- 

Since all the Bethe roots such that l^{Aa) € H 1; ]] are contained in where ^ is increasing and owing 

to minima) < n and maxj^Q) > - n, it follows that 

[f;^lc?(A). (3.21, 

Moreover, for L large enough 

?(/a) c {x € R : d{x, ^(/a)) < II?- c ^(/a+./ 2 ) . (3.22) 

Then, passing to the limit in (13.201) . it holds that 

[0;D] c^(/a+,/ 2 ) ^ [-v;D + v] c^(/a+,) (3.23) 

for some v > 0 small enough. This inclusion ensures that 

c ■ (3-24) 

Thus, due to (13.201) . the points 

?L = and qR = (3.25) 

are well defined and belong to 5;^ ^(/a). Moreover, when A > 1, one has the bound 


qR-qL<^ (3.26) 

as ensured by the fact that ^ is strictly increasing on I^+e and satisfies ^{x + n) - ^(x) = (L- N) IL < NIL owing 
to 2N < L. 

Furthermore, since ^{Srj^^ilA)) is an open L-independent neighbourhood of the compact interval 

[-v/2;D + v/2]d[^; ^^^^^^ ], (3.27) 

by compactness, there exists a > 0 such that the domain G{a) as depicted in Fig.[T]satisfies G{a) c ^(iS;,,e(/A)). 

Finally, for any ajL e ^(•S;;^£(/a)) it is possible to define Xq = ^“^(a/L). This holds, in particular, for a = hp, 
p e H 1; n ]]. Furthermore, owing to |/la| < A, if £ [[ 1; 77]] the strict increase of ^ on I a ensures that 

{Aa : 4 e 11; iV]]) - {{x,}f \ ■ (3-28) 


Finally, define = {Aa}^ \ {{xq}^ \ {x/j^)j|. One needs to recourse to such a definition since it could be that 

PalL t ^(5p,e(/A)). However, for those pajL e ^(^^,^(4)) one does have Xp^ = ^~\palL). 

These properties being established, it follows from a straightforward computation of residues and the strict 
increase of ^ on /a+e that 


^(4) = 


P(4) 

2n 


2nL 


2 [0(4 - XpJ - 9{A - Xftj] + ^ 


9(A-pX'(p) dp 

Q2mI4(R) _ I ItT 


(3.29) 
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where the contour C is defined as C = ^ U r_) and r+ have been depicted in Fig. [U The expression can be 
further rearranged. Namely, setting Ce - one has 




c 


2n 


m 

-J + 

?L 



6(A-/u)e(M) 

Q-linemii) _ \ 


In 


m 

- -^[e{A-qR) + e{A-qL)] - JK(A - ju)?sym(M) ■ dju + ri[?](T) . (3.30) 

?L 


In the second line, appears the function 

4^sym(d) = ^(A) - 

and ri is the remainder 


- - Z f 

e=±J“ 

r. 


K(A-r\z)) 


• In [l - ■ 


dz 

2inL 


(3.31) 


(3.32) 


Note that the second line is obtained by carrying out an integration by parts followed by a change of variables in 
what concerns the form of the remainder ri [^](/l) given in (13.1071) . 

Thus, all-in-all, one gets the representation : 


qR 

?sym(d) + f^(d- M)Isym(M) ' ^ ^ -?l)) 

?L 

1 " — 

- ■ (3-33) 

a=l 

The latter is already enough so as to characterise the limit Indeed, define qi = ^“^(0) and qR = Then, 

one has 


\qL-qL\ < + |r'(:^)-r'(o)| 


C2, 


"2L' ’ "2L'' "2L' 

< 0 ( 1 ) . ( 334 ) 

Nofe fhaf, in fhe second chain of bounds, I used fhaf 


iir-r 


llL‘”(f(5,,,(/A))) 






(3.35) 


as ensured by Proposition lA.ll Similarly fo (13.341) . one concludes fhaf \qR - qR\ = o(l). These estimates broughf 
fogefher wifh (13.261) ensure fhaf, for A > qR - qL < z:. 

If remains fo bound ri[^](/l). For any d G 1. one has 


ri[fl(4)| < !«..(*,«) ■ I, If'-g/ 


1^ < r 

2nL ~ L2 ■ 


(3.36) 
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The last bound follows from inf 52 , 2 e(tA) 1^1 > \ > 0 for L large enough. Also the XjlJ decay in 

(13.361) follows after an asymptotic estimation of the integral by a variant of Watson’s lemma : the boundaries of T^ 
generate an algebraic decay in L starting with 0(L“^) while all other contributions to the integral are exponentially 
small. The uniform convergence of to ^ - Djl and the above bounds are enough so as to take the 

L —> +00 on the level of (13.331) . This yields the system of equation for three unknowns : the function ^sym 

(id + ;,,])[^sym](d) = ^ - - qR) + - q^)] (3.37) 

and the endpoints of integration qi, q^ 

^sym(^s) = ~^sym(^L) — • (3.38) 

Also, the additional condition ^ ^ is imposed for A > 1. 

In virtue of Proposition 11.11 there exists a unique solution to (I3.37I) - (I3.38I) given by qn - -qi = q and 
^sym(d) = p{A I q), with q the magnetic Fermi boundary associated with D. This ensures the bounds (13.141) and, 
upon elementary manipulations, the form taken by the non-linear integral equation (13.151) satisfied by 


• Diffeomorphism property 

Lemma fTTSl ensures that ; D]) = [-^;^. Furthermore, since D ^ D with either D < 1/2 for -1 < A < 1 
or D € [0; 1/2] for A > 1, one has for L large-enough that 'q <2q < -i-oo. Then, by Proposition lA.ll one readily 
deduces the biholomorphism property from (13.141) . 

Then straightforward bounds in the non-linear term based on the fact that 


sup {|5i/?/^(/l,/r)|-H|/?/^(/l,/r)|} < C-g^{A) with 


( 1 

for 

A > 1 


gA(4) = 1/(42+ 1) 

for 

A - 1 

(3.39) 

( 1/cosh(24) 

for 

-1 < A < 1 



ensure that 

?(T) = ^'(T|?) + 0(^gA(4)) 


(3.40) 


with a remainder uniform in L and 4 e R. Since ^a/^q bounded on R, this ensures that, for L large enough 
> 0. Therefore, ^ is n strictly increasing diffeomorphism from R onto ^(R). The explicit form of the range for 
-1 < A < 1 follows from computing the limits of ^(4) at 4 ^ ±oo starting from the definition (13.11) of the counting 
function. When A > 1, the finite difference growth ^{x + n)- ^(x) = (L-N) IL ensures that ^ is a strictly increasing 
diffeomorphism form R onto R. 


• Uniqueness of solutions to the non-linear integral equation 

By hypothesis, for -1 < A < 1, there exists 0 < v small enough such that l/2-D-v>0 and 

^ (^-£>-v)); ^(^-D-v) + d] for a = l,...,n. (3.41) 

7T ^ JX ^ 

Then, by Lemma [TTSl for -1 < A < 1 one has \'q) a ly for some y > 0 and uniformly in L. Likewise, 

when A > 1, given C > 0 such that |pa/L| < C, one has \ 'q) ^ ly for some y > 0 and uniformly in L. 
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Suppose that one is give two solutions ^ 1,^2 to the non-linear integral equation (13.151) and satisfying to all 
the other requirements, (13.141) in particular. In virtue of Proposition lA.ll there exists 77 , e > 0 such that for any 
z e I ?), one has 






A-CW dd 

and 


r(4) - Z 


_ 7 2i7r 


inf inf \^a{A) - z| > cq > 0 

AedS2r,,2AIy)zeSr,AIy) 


(3.42) 


dS2q,2€Uy) 

for a = 1,2 and some constant cq only depending on ^o(* I <?)■ Likewise, one has 

II?;' - I ^\\L«{Us,My)\‘i) ) ^ • l?« - I ?)IIl~(5.,(/27)) « = 1’ 2 , 

and, after straightforward bounds in (13.421) . one gets that, for some constant C > 0, 

11^1 -^2 llL“(fo(5,,,(/7)|^)) - ‘^•ll^l ~ ^^h‘>^{S^Jl2y)) ■ 

I stress that, in virtue of the condition satisfied by the p^’s, provided that L is large enough, for all 

11 { ^ 

|/ji,...,/j„,pi,...,p„}u [ 0 ;A-f 1 ] itholds - £ ^o(*5;,,e(/y) | ^) . 


(3.43) 


(3.44) 


(3.45) 


All is now in place to estimate the norm ||^i - ^ 2 11 


(/y)) 


by using the non-linear integral equation satisfied 


by ^a- Agreeing upon X{-^a = in virtue of (13.441) . it holds 

Cn 


4 (* I i'^P.An ; 4 (* I ; {Xh^^n) 


L-^iS^Ahy)) 


11^2^^11^00(5^^ (/2y)XR) ■ 11^1 ^^h<^{S,^{l2y)) • *^ 2 - 46 ) 


Further, define 

?R;« - and qL,a - ■ 


In order to bound 91v;2[ ^](d), it is enough to observe that, in virtue of (13.431) . 


C 


< Y and ?R;l-^i?; 2 | < C ■ \\^\ 
Then, it holds, for any A e Sx/^ihy) 


_ 9 ^' 

HS^^ihy)) - L ' 


(3.47) 


(3.48) 


9iv;2[^l](d) - 9lAf;2[^2](d) 


‘iR.i 9K;2 

J Rj/A,m)IIi(m) -? 20 u)]dp + J Rj/A,ju)la(M) -? 2 (?i?; 2 )]dp 


9K;1 


< 


ll^^?llL”(52,(Py))-{^-^«4Hl^l-^2lli.(5^^(,^^)) + f?fi;2-^«;lfll^2llL~(5,,(/y))} ^ - ^A-{s.^il2y)) ' 


An identical type of bound can be obtained for 91v;3- Finally, using that, for L large enough and independent of a 




>0 


(3.49) 
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one bounds 3?^^;! as 



(3.50) 


By taking the difference of the two non-linear integral equations satisfied by and ^ 2 ; the various bounds obtained 
earlier lead to 



(3.51) 


The latter can only hold for L large enough provided that thus entailing uniqueness of solutions. 

• is not necessarily positive on [-A; A]. 

In this last part of the proof, I study the case where, a priori the limit of the extracted subsequence is not 
positive and show that such a situation cannot arise. 

Prior to going into the details of the analysis, one should observe that this situation cannot arise for -1 < A < 0 
since then one has the trivial bound 2n^'{A) > p'(/l) > 0 on /a. 

When A > 1, it is readily seen that the functions and its limit are n periodic. In principle, one could have 
that A > njl. Yet, then it is enough to observe that the Bethe roots can be decomposed as Aa - Aa + ria n, where 
Ha £ Z is such that Aa £] - n'/2 ;7r/2]. Making explicit the dependence of the counting functions on the Bethe 
roots, it holds 


I {d,}f) - ^'{A I {Tjf) . 


(3.52) 


Therefore, since both functions will admit the same limit of the extracted sequence, it is enough to reason on 
the level of the Aa which all belong to the interval [-7r/2 ;7r/2]. Once that it will be established that > 0 on 
[-7r/2 ;7r/2], its strict positivity on intervals of large diameter will follow. Hence, below, when A > 1, we shall 
assume that 0 < A < njl. 

The function is holomorphic on Sxf,{h\), and thus admits a hnite number of zeroes on [-A; A]. Let 




(3.53) 


be such that 


[z £ [-A; A] : ^\z) = 0 or z = ±a) = 


(3.54) 


Also, given d > 0 and small enough, let 



(3.55) 


By Proposition I A. 1 1 there exist e,T] > 0 such that 


^ : ‘S2,,2.(/f) ^ ^{32,MI?)) 


(3.56) 


is a biholomorphism that satishes 




(3.57) 
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and 




(3.58) 


is a biholomorphism. Furthermore, one has sgn(^') = on 

Let m® he the multiplicity of the zero 3 ®. By Rouche’s theorem, since \\^' - ^ 0, given any 

d > 0, there exists Lg such that for any L > the function has zeroes, counted with multiplicities, inside 
of the disk Dj® ^^ 4 . The distinct zeroes will he denoted hy with a - 1,... A priori some of these zeroes 
can have non-zero imaginary parts while other will he real. Thus, I assume that the zeroes are ordered in such a 
way that 

for a=74^^ + l,...,7i(^) . 




and 


Further, take L large enough so that 

[1(3® + d) - - d)| > ^1^(3® + d) - ^(3^*^^'^ - d)| > 0 . 

Then, define ^®^ e 7® as the two solutions to = -1 such that 


_ ^{k+i) maximal 


and 


3^^^ -I- q''j^’ is minimal . 


rik) 


(3.59) 


(3.60) 


(3.61) 


In virtue of (13.601) . for L large enough, both < 7 ®^ are well-defined, exisfs, and are distincf. Furthermore, for any 
given d > 0, there exists L large enough such that 


k=0 

Finally, let 


ZIK*’- 




-I- 


|?l^^-3® 


< 4(r -t l)d . 


X = (x£[-A;A] : - l) , = Xn{|j/®) , - X \ X' 


(in) 


k=o 


Similarly, setting Y = (Ta)^, 

yhn) ^ y J5f(in) y(out) ^ y \ T^*") . 


(3.62) 


(3.63) 


(3.64) 


The purpose of the below paragraph is to estimate The first step consists in estimating the number of 

roots in the interval with k = - 1 ,... r under the convention that 


■ 5 ^ = -A and = A . 

has constant sign on each of the intervals 

”Re 

Due to the ordering (13.591) of the real roots one gets the upper bound 


(3.65) 


(3.66) 




+ 
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In the first line, [*] denotes the integer part. Furthermore, so as to get the last bound, L is taken large enough so 
that ^ Then, summing over k, and using (13.621) . one arrives to 

< 8(r+l)dL-|ri|^„(^^) . (3.68) 

In order to write down a non-linear integral equation satisfied by I sfill need fo define auxiliary contours. 
Lef G®(a) and ifs boundary u be defined as in Fig. |2] By compacfness, if follows fhaf fhere exisfs a > 0 
such fhaf G^^\a) c for any k = 0,... ,r. 



Figure 2 - Domain G ^^\a) delimifed by ifs boundary F® UF^^ When - 1 one has p ® and ^ 

while, when = -1 one has and ■ 


If now remains to characferise fhe sum over Befhe roofs occurring in The laffer can be recasf as 


-^^K{A-Aa) = 2 ] K{A-y) " K{A-y) 

a=\ yEF<™') ygyCin) 

^ V m I J. ^ V m 'I V ~ 


ygyCout) 

where C® = U F®). Now, agreeing upon 

J\p)KiA-p) 




e2-^(/')-l 
c® 


dp (3.69) 


C+^ - n H-, one can decompose the integral over C® as 

(ju)K(A-p) 


—(*) 


■dp = -/r® r KiA-p)^'(p)-dp + /r® V e f — 

Qlmmp) _ 1 J ^ j - 1 

c® 

,(*+ 1 ) 


dp 


= r K(A-M)^'(M)-dp + r®[^](d). (3.70) 

VO 
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There, I have introduced 


‘Ir 


■(« 




dz 


(3.71) 


‘?L 


W J»+ 1 ) 


'(« 


Finally, introduce the intervals 

j(±) = ^ 

k : /f®=±l 


(3.72) 


All the above leads to the non-linear integral equation satisfied by ^ 

(id K/+) -Ky(-))[^'](/l) = - (poutK'^J + (PmK.yi) + 

^ k=Q 

There, I have set 


- <^out(d) + 0i„(d) + j]r®[^](d) . 


0out(d) = \ Y, K{X-y) and 0in(T) = j ^ K{A-x). 


(3.73) 


(3.74) 


ygJ^(OUt) 


JC€Z(i")\y(m) 


Let R/+) be the resolvent operator@to id-i-Ky(+). By Lemma fOl the resolvent kernel is strictly positive : Rj(+){A,n) 
0. Furthermore, by (11.161) 


> 


P;„(d) = (id - R^„)[^](d) > (id - LR)[^](d) - poo(d) >0 on R . 
Thus, adopting the notations of Lemma [L4l one gets 


(id -[R/o]/-,)[r]W = Pji.M) - 'Aout(d) + .Ain(d) + j]r®[^](T) 


k=0 


where 


<Aout(d) ^ ^ ^ Rj(+M,y) and i/^in(T) = ^ Y Rjm{A,x) . 

ygy(out) xeX(‘")\y<‘") 


(3.75) 


(3.76) 


(3.77) 


Finally, 






J 


7?/+)(d,ju)^'0u)-dyu -t 


^ — 4 - 


Rj(.M,r\z)) 

QlmeLz _ I 


• dz 


(3.78) 


■z^k) ,(*+ 1 ) 

It 


~{k) 


It follows from Lemma [L4l that the operator id - [R/c+j]/-) is invertible and that its resolvent operator [.^^/+)]/-) 
has a strictly positive integral kernel [3ij(+)\j{-){A,^) > 0. Therefore, one has 


r(T) - (id + [^/o]^(-))[p;„ - <Aout + <Am + j]rd^)[^]](d). 


(3.79) 


A:=0 


f. Recall that when A > 1 it holds diam(y*^' U 7* 9 < so that indeed the involved operators are invertible in virtue of Lemma [T3] 
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By the above, the operators and afortioti id + Mj(-) are strictly positive. Since, > 0, lAinCd) > 0, this 

leads to the lower bound 


r(d) > P'jM + (id +[^/w]y(-,)[ - ^out + 


(3.80) 


k=0 


Since y ^ and pi yf(out)j bounded by the bound (13.681) ensures that 

_ ^y(out) ^ _ 

lAout(d) < II-^j(+)IIz.o<.(r 2 ) 2 j^i - leadingto limsup |(id+[^y(o]/-))[ - iAout](d)| <C'S (3.81) 




by continuity of (id + [^/(+)]/-)) and with a bound that is uniform on R. 

Owing to ||^||l“(/a) ^ 'A\^\\Lr(itd ^’^dG [-A; A], it follows readily that the contours F® are all bounded 
in L and thus for any A e I\ 


i :=0 e=± - 


f 

'm 


/?j(o(d,rAz)) 

g2iffeLz _ 


■ dz 


< 


WRji^A 


‘’(/aX5; 


"a'- 


c 

L 


(3.82) 


where the 1/L decay follows from a straightforward estimation of the integral while the bound on the integral 
kernel Rji+) follows from the inclusion c what leads to that ^“^(rji) c S 2 -q, 2 e{lf^) c S,if^{I\+ 2 e) 

as ensured by (13.571) . Thus, in virtue of (13.621) . one gets 

r 

limsup (id + [^/o]/-) - 0(d). (3.83) 

L^+oo 

This being settled, it remains to take the L —> +oo limit superior of the inequality (13.801) followed by sending 
d ^ O'*'. One gets, for any A e I\ 

f (T) > P'M ^ Poo(d) > inf Poo > 0 . (3.84) 

This contradicts that, either ^'(A) < 0 on /a or that has zeroes on /a. ■ 


3.2 Various corollaries of interest 

Theorem 13.21 has several important corollaries. First of all, it guarantees the uniqueness of solutions to the 
logarithmic Bethe Ansatz equations (10.71) throughout the regime A > -1 and, in particular, for A > 0 when 
convexity arguments cannot be used. To the best of my knowledge, uniqueness of such solutions has never been 
proven earlier. 

Corollary 3.3. Under the hypothesis of Theorem 13.21 and point v) in particular, there exists Lq such that the 
system of logarithmic Bethe Ansatz equations associated with the given choice of integers {£a}^ admits a unique 
real-valued solution for any L > Lq. 

Proof — 

If there would be two solutions to the logarithmic Bethe equations associated with the given choice of integers 
than one would then be able to build two distinct counting functions a = 1,2. These, however, will 
both satisfy the non-linear integral equation when L will be large enough. As follows from Theorem [T2] point v), 
the non-linear integral equation admits a unique solution, contradicting that^^^ for L large enough. ■ 

A second important consequence of Theorem l3.2l is the existence of an all order large-L asymptotic expansion 
of the counting function. 
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Proposition 3.4. Let hi < • • • < /j„, e |[ 1; ]] and p\ < ■ ■■ < Pn , Pa € Z \ [[ 1; ]] Ae and increasing 

sequence of integers such that 
• for -1<A<1, A = cos(^), 


n-q i _ N'-\ . pn-N' Pi-I 7T-q\ _ N'-\ . 
n ^2 L L ’ L ^ n L >' 


(3.85) 


• for A > 1, the integers {pa}^ are uniformly bounded in L : \palL\ < C for some L-independent C>0. 

Let be the solution to the logarithmic Bethe equations subordinate to this choice of integers and let 

N' - N + s with s bounded uniformly in L. Then, the associated counting function ^ admits the large-L asymptotic 
expansion 


f{A) - m\q) + 

p=\ 


(3.86) 


There Xa are defined through f(xa) = alL. 

This asymptotic expansion is such that the remainder just as all functions, ^p, p - 0,... ,r are holomorphic 
on iS^/ 2 (K.). Furthermore, the remainder is uniform on S^for any 0 < q < ^jl. The function A as defined in 
(11.511) with'q the solution to the magnetic Fermi boundary problem associated with D = NIL. Also, 


= ^~\A\{Xp’J[*'-,{XhXl) [Pa]T' = [PaTl 


- roA)/ 


U|aK, 


(3.87) 


given in (13.171) and 


ff{A\{xpJf,{%J\) = 


1 


2 ^n(?l q) 


eRi^iA, eq) ■ (^f{eq \ {Xpjf, {xhXi) ~ ~ ^ 


(3.88) 


Finally, the unique solutions qi/p to f(qL) = 1/(2A) and ffqp) = (N l/2)/L admit the large-L asymptotic 
expansion 


r Jk) 

— — V 

qR-q - 2^-jT 

k=\ 



and 


qi + q 


r 






(3.89) 


where the first terms of the expansion take the form 

ll2-^^\±q) 


ql^ = 


q) 


(3.90) 


Also, the second order deviations in respect to q read 
1 


( 2 ) 

q± = - 


{m\^Y 


^(*9 I 9) 


1/2-4^\±?)] (^;^^)'(±?) + ^f{±q)f'Q{q I ?) + [1/2-^«(±?)]2} . (3.91) 


The dependence of and on the auxiliary rapidities has been omitted both in (13.901) and (13.911) . 

The decomposition N' = N -v s might appear slightly artificial in that one can absorb all the dependence on s 
by simply doing the substitution {D,'q, s) {D','q',0) where D' = N'IL and^' is the magnetic Fermi boundary 
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associated with D'. However presenting the expansion in the above form allows one to immediately take into 
account various terms of the asymptotic expansion when D —> D sufficiently fast in L for instance 

D-D = o(^) while D'-D = o(i). (3.92) 

Further, note that all the terms of the asymptotic expansion depend implicitly on 'q. It is readily seen that the 
latter admits an all order asymptotic expansion in D - D. Presenting the asymptotic expansion of ^ in the form 
(13.971) allows one to absorb all the orders in D-D and thus have a slightly simpler to obtain, from the computational 
point of view, asymptotic expansion. If one would make the additional hypothesis D - D = then one 

would also have'5^- q = and it would be possible to simplify (13.971) further by replacing '^^by q. The 

presence of terms D - D in the asymptotic expansion was first noticed in If34l . 

Proof — The first two terms of the asymptotic expansion, just as its form, are readily obtained by recasting the 
non-linear integral equation (13.151) in the form 

?(d) - ?o(^ I ?) + I r'; + 9?;v[?](d) (3.93) 

and then using straightforward bounds on the "remainder" operator 9 ?a?. The starting point for pushing the asymp¬ 
totic expansion one order further, viz- up to r = 2 , is to observe that the properties of the solution to the non-linear 
integral equation lead to an overall estimate 

[?-?«! + |?+?l| ^ 0(i/l) . (3.94) 

A straightforward application of Watson’s lemma to 9?]V;i[^] yields 


9?W;l[^](4) 


1 (7?/^(d,^/?) 

4;r2L2l 


?(?l) 




1 

2AL%fq\q) 


(Rl^{A,q) 


-Rl^{X,-q)) 



Also, due to (13.941 ). it holds that 

9?yv;2[?](d) + 9iiV;3[?](4) = 7?/^(T,?p)?(?p)^ 

The above bounds already ensure that ^ admits the large-L asymptotic expansion up to 0{Lr^) : 


ad) = a(d|?) + {xpX,-{xhJ\) + 0{L-^) 


(3.95) 


(3.96) 


(3.97) 


where a and are as defined in the statement of the proposition whereas the expression for involves the 
L-dependent endpoints 'qRjL '■ 


6 =± 


Ri^{A, eq ) 
^24^'(^|?) 


2 


[Ri^{A,q)(qR -'qf 


Ri^{A,-q)(q + qif] . (3.98) 


To conclude, it solely remains to obtain the asymptotic expansion of < 7 ^ - ^ and qi + q vrp to 0(L 2). The latter 
can be obtained by expanding the defining relation 


1 -v-x A^-tl/2 

— - a^z.) and - - — - ^{qR) 


(3.99) 
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to the second order m'qR-'q and 1^1+ 'q and using the a priori estimates (13.941) so as to separate slower and faster 
decaying terms. All-in-all one obtains the expansion (13.891) up to the second order, viz. r = 2, where is as 
defined in (13.901) . However, at this stage of the analysis, the expression for < 7 +^ has still not its final form in that 
it is given in terms of To conclude, one should first insert the expansion (13.891 ) to the first order in 1/L into 
(13.981) so as to get (13.971) with r = 2 along with the form of the coefficients (I3.87I) - (I3.88I) . Starting from there, one 
readily obtains the claimed form of the second order coefficient ^ as given by (13.911 ). 

The existence of the all order asymptotic expansion is obtained by a classical bootstrap argument, which I 
shall not detail here. I refer to Section 3.2 of l28l for similar handlings. ■ 

3.3 The counting function when D = 1/2 and -1 < A < 1 

Theorem 3.5. Let —1 < A < 1 and N/L - D - 1/2. Then, for any e > 0 there exists Lq such that the counting 
function f constructed from the Bethe roots associated with the choice of integers €a= ci in the logarithmic Bethe 
equations admits the expansion 

A 

^(T) - ^o(d I -^oo) -I- 0(e) with ^o{A \ - 1 - 00 ) - J'pooip) ■ dp (3.100) 

0 

where the remainder in uniform in L > Lq. 

Proof — 

Let {/Iq) be a solution to the logarithmic Bethe equations associated with the choice of integers 7 ’q = a. In 
virtue of Proposition 12.31 there exists A^ such that < e, where cX has been defined in (12.561) . 

The counting function f associated with these roots is a sequence of holomorphic function on .S;^^(R) and, as 
such, admits a converging subsequence to some f. It remains to characterise the limit f and hence the limit of ^ 
according to the reasoning of the proof of Theorem l3.2l 

The first stage of the proof follows the analysis developed in the proof of Theorem [3]2] relative to waiving-off 
the possibility that f has several zeroes on R. Taking the same definition of (13.631) with A now being replaced by 
Af one arrives to (13.791) with iXin/out defined exactly as in (13.771) . Again, by positivity of the integral operator and 
by invoking (13.751) . one gets the lower bound (13.801) . viz. 

r 

?(d) > Peo(T) + (id +[^/o]y(-))[ - ^?out + 2]7®[?]](d) . (3.101) 

k=0 

By continuity of id - 1 - [l%j(+)]j{-), it follows from < e and the bounds (13.681) on which allow one to 

control the cardinality of p| yf(out) 

fimsup|(id +[^jw]ji-))[ - iAout](d)| < C • {d(r-i-l)||f -h e} . (3.102) 

L—»+oo 

In the above formula, r refers to the number of zeroes of on ] - A^ ; A^ [ and thus depends a priori on e. One 
also obtains a similar bound on the second term, namely 

r 

limsup |(id + [^y(o]y(-,)[ V7®[?]](d)| < C' • d • (r + 1 ) • |irilL“(/^p - (3.103) 

L^ + 00 

Taking the pointwise in T L —> - 1-00 limit superior of (13.1011) . one finds the lower bound 

f\A) > poo(d) - C"{6 • (r + 1) • |irilL~(/^p + e) . (3.104) 
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It then remains to send first d —> 0^ and then e ^ 0“^ so as to get the strict positivity of ^ on K.. 

Now, taking for granted that ^' > 0 on R, one again picks e > 0 and such that ca^ < e and definesresp. 
fo be the closed to A^, resp. -A^, solution to exp [2mL^{A)] - -1 lying outside of /a^. Repeating similar 
handlings to the ones carried out in the corresponding section of the proof of Theorem [321 one gets 

?sym(d) + J K(A-jU)Isym(M)-dM - ^ 

?L ;6 

- ^(?sym(?R;.) • 0{^-qR-,e) - ?sym(?L;.) ’ -?L;.)) + ri[?](d) - ^ 9{A - x) . (3.105) 

x€y(™» 

Here = {Aa : Aa i ;?R;i 5 ]) and ri[^](/l) is as given by (13.1071) relatively to the above defined endpoinfs 
of infegrafion^R/R;f. Affer a few handlings, one recasfs fhe above equafion as 


^symCd) = p{A I +CX3) + Z{A I +oo) 


2n 


{^sym(5i?;e) ^sym(5i;e)} 


+ - 7 2] ip{A-x) (3.106) 


^gy(out) 


where 


9? 


fj^[^](T) = 9?A;l[^](d) + J R(A-ju)l^, 


(tot) 

N: 


symifO 4'sym(<?R;6)] ' d/T 




J R(A-M)ll 


symiu) ^sym(?L;.)] • d// (3.107) 


and 




^2meLz 


dz 

2inL 


(3.108) 


If remains fo bound the various terms in (13.1061) . A counting of solution argument ensures that 

[L(^sym(-Af) - ^sym( “ “))] < and [^(^symC + 0°) “ ^sym(Af ))] < Lca, ■ (3.109) 


Thus, 


^sym(5i;e) ^sym( ‘^) ^ fsym( A^) ^sym( '^) ^ ^A^ + 1/T 


(3.110) 


and a similar bound holds for ^sym(+ °°) - ^sym(9R;e)- From there, since ^sym(+°°) + ^sym(-“) - 0, one readily 
bounds the second terms in (13.1061) . From the above bounds one also infers that 


+ 0O 

J R(A-M)ll. 


sym(p) ^sym(<?R;6)] ■ d/r < ll-^llLhR) ’ (+ 1/F) 


(3.111) 
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and analogously for its counterpart related to Finally, one gets that 


9?W;l[f](d) 



(2Ae + 1) • 

inf{f(5) : 5e5;,(/Aj} 


(3.112) 


Given e > 0 fixed, the infimum is bounded from below due to the convergence of f on compact subsets. This 
entails the claim, upon taking L large-enough. ■ 


4 Applications 

In this section I establish several more or less direct applications of the existence of the large-L asymptotic 
expansion of the counting function. The first of these corresponds to the proof of the existence of limits of the 
type (10.131) . The second application concerns the proof of the conformal behaviour of the spectrum of low-lying 
excitations above the ground state in the so-called massless regime of the model, namely when there exist zero 
energy excitations. This regime corresponds basically to D €]0; 1 /2[ for any A > -1 and D = 1/2 for -1 < A < 1. 
Below, I shall confine myself fo fhe regime 0 < D < 1/2 since fhe laffer is fechnically much simpler fo deal wifh. 


4.1 Densification of Bethe roots 

Theorem 4.1. Let N,L -i-oo in such a way that N/L — > D e]0; 1/2]. Let q be the unique solution to (11.251) 
subordinate to D and let the particle-hole integers (1231) satisfy to the hypothesis of Theorem \3.2\ Then, given any 
bounded Lipschitz function / on R with Lipschitz constant Lip[/], it holds 

L 

\ ^ f I (4.1) 

«=i 4 

where {Taj^ correspond to the unique solution to the logarithmic Bethe equations subordinates to the given choice 
of particle-hole integers. 


Proof — 

Due fo fhe unbounded nafure of fhe Befhe roofs at D - Ijl and -1 < A < 1,1 firsf deal wifh fhe simpler case 
of a uniformly bounded in L disfribufion, namely when 0 < D < 1/2 and -1 <A< lorO<D< 1/2 when A > 1. 

For fhis range of D and A, Proposition |23] ensures fhaf fhe Befhe roofs such fhaf ^(2^) = a/L all belong fo fhe 
compacf [-A; A]. Also, if has already been esfablished fhaf, provided L is large enough 


for some e,q > 0 small enough. Finally, for L large enough, if holds o/L e ^o(^a+6 I <?) for any a e [[ 1; fVI- The 
finile sum one sfarfs wifh can be decomposed as 


1 ^ 

t 2 ]/(^) = + ^2 + ^3 + ^4 


wifh 


•S'! = 


a=l 


0=1 


while 




(4.4) 


0=1 


0=1 
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and 


^4 = 

a=l 

, 5^4 is a Riemann sum and, as such, converges 
D q 


(4.5) 


=^4 


J/ ° ^0 ^(‘^ I <?) • Jf{X)p{^ I ^) • dd 


(4.6) 


where the last equality follows from ^o{\_-q ',q\ | < 7 ) = [0; D] and ^'^(A \ q) = p{A \ q). Further, Rj^(A,p) is readily 
seen to be smooth in or and satisfy < C for a bounded. The latter leads to 


^3 < C' • Lip[/] -Iq-ql-y 0 . 

L, L—>+oo 

Clearly, ^2 = 0(L“^) since / is bounded while 


N 


hi < j-Lip[/].|ir‘-f„-‘(.i?)ihfe(,^„i,)) = o(L->). 


(4.7) 


(4.8) 


I now focus on the remaining cases, viz- D = Ijl and -1 < A < 1 and picking some M > 0 decompose the 
sum of interest as 

7 2/0.) = 'j 2 /MO + 7 2 [/“/'‘(t) - I +“)] + 7 2 /°/o’'(7 ' +“) ^ 


a=l \Aa\>M a:|4a|<M 

The last term converges as a Riemann sum 

M 

mpum^ • 


|4J<M 


I Z ^ /- 


(4.10) 


The estimates of Theorem [33] ensure that, provided e is small enough, one has 
limsup i I+oo)]| < C-Lip[/]-e. 


(4.11) 




l^|<M 


Finally, 


limsup j V f(Aa) < II/IIl-’CR) • limsup(|#|a G E 1; : |do| > M]) . 

T . __ /, • * r . I_, d 


(4.12) 


L—>+00 


\Aa\>M 


L —^+00 


Thus, upon sending L —> +00 and then e —> 0'*' one gets 


M 


^ J /(d)poo(d)dd + 0 (limsup(Z{a€ll;iVI : |TJ > M})) . 

-M 

Finally, it remains to relax M — > +cx). The remainder term, in virtue of Proposition 12.31 tends to 0. 


(4.13) 
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4.2 The conformal structure of the low-lying excitations 

The XXZ chain embedded in an external magnetic field h refers to the Hamiltonian Ha;/i = Ha - on hxxz- 
Its eigenvalues and eigenvectors are readily deduced from those of Ha owing to [Ha, S^] = 0. The eigenvalue of 
Ha;/j associated with the eigenvector 'P({/la}^) of Ha that is parametrised by the Bethe roots {Aa\^ takes the from 

h ^ 

£({T,}f) - (7A - -)l + 2]e(T,) (4.14) 

a=\ 

with e as defined in (I1.49I ). The large-L expansion of the above eigenvalue involves the effective dressed energy 
defined as 

eeff(4 I Q) = e{A I Q) + s{Q I Q)(^{Q, A\ Q) A \ Q)) . (4.15) 

Proposition 4.2. Under the assumptions ofTheorem \3.2\ for any 0 < D < Ijl and A > -1 /t holds 

2 . 

= L ■ £o(?) + 2] ixp'jT ^; + 0{l-^ ) (4.16) 

k=i ^ 

where 


£o(?) ^ \-(p)P(h I?) -dju (4.17) 

«+5 n 

£l(?l ^ I?) - \q) (4.18) 

a=\ a=\ 

with fieff tis defined in (14.151) and, finally, 

■Axhfi'l) = + (4-i9) 

with as defined by (13.901) . 

Proof — 

Since e is holomorphic in a neighbourhood of (he real axis, one has the representation 


h r ” 

e({4fl}f) - (/A - -)l ^ L I tiifAip) -dp + ^ [e(ApJ - ^(xtifi] 

- a=l 

qt 


In [1 - ■ — . 

2i7r 


_ y f i„ r. _^ 2 meLz- dz 


A straightforward application of Watson’s lemma leads to 

Ar, d(q) 

in 11 — e I • t:— =- 

f'if-Kz)) 


y f In [1 - ■ — = 


2 i 7 r \2Lf'(p\q) 


+ 




(4.20) 
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In its turn, the one-fold integral admits the expansion 


J e(ju)f(yu) ■ dyU 

?/. 


2 l/i 

+ o(^) 

/t=o ^ ^ 

qL 


(4.21) 


where, so as to lighten the notation, I have dropped the dependence of the ^^(4) on the roots Xp^/h^. It then remains 
to use the Taylor expansion 



(^+ 1 )! 



(4.22) 


and the differential identities satisfied by the dressed phase (11.531) so as to get the claim after some straightforward 
algebra. ■ 

Lemma 4.1. Let he > h > Ofor -1 < A < 1 and he > h > h^p for 1 < A, then the function Q i—> Sq{Q) attains a 
unique minimum at gp £ 4 the so-called Fermi boundary of the model. Furthermore, Qy is such that 


£(4 I 2f)|/2p <0 , e(4 I 2f)|/^p > 0 and £(2f I 2f) ^ 0 . 


(4.23) 


Proof — 

It follows from straightforward algebra that 

(&T){Q) - £{Q\Q)-p{Q\Q)- (4.24) 

Since, p{Q\Q) > 0, the derivative has the sign of e{Q \ Q). It follows form Proposition 1 1.21 that s{Q | 2) < 0 on 
[0; 2 f[ and e{Q | 2) > 0 on This entails that Q i-> £o(2) admits a unique minimum at 2 f. ■ 

Define the Fermi density Dp by Dp - p{Qp \ Qp). Given L large enough it seems reasonable to expect that 
the ground state of with N such that \Dp - NIL\ is minimal will give rise to a0 ground state of In the 
present state of the art, I am however not in position to prove the statement due to the lack of a sufficient uniform 
in L, N control close toNjL - 1/2 of £({4a}^). 

Theorem 4.3. Let h be such that Dp £]0 ; l/2[ and assume that N, L are such that Dp — NIL = 0{L~^). Let s be 
fixed and consider the solution to the logarithmic Bethe equations subordinate to the choice of particle- 

hole integers {pa}1 and {ha}" 


Pa = 1- Pa for a=l,...,np and p„-^+a = N sPa 

for <3=1,. 

..,«p , 

(4.25) 

with Pa € N* fixed in N, L, 




- h~ for a = 1,... ,nj^ and h,q+a = N -t- s -t- 1 - ha 

for a - 1 ,.. 

■’K ’ 

(4.26) 


with fixed in N, L, and n^ -i- n^ = n = n* + n^. 

Then, it holds 


f {- 4 + I Qri -1] + 

"p "I 

+ + y,\-k-v2]) 

e=± a=l a=l 

t. There could, in principle, more than one ground state since \Df — N/L\ could attain two minima for some fixed value of L. 


) + o(^) 
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where 


^ s\Qf I Qf) 

P'{Qf I Qf) 

refers to the velocity of the excitation lying on Fermi boundary. 


(4.28) 


Proof — 

Observe that the choice of integers given in the statement of the theorem always gives rise to a solution of the 
logarithmic Bethe equations since the conditions stated in Proposition 12.11 are always verified provided that L is 
large enough. Further, since D q{D) is smooth, it holds that Qf - ^ = 0(L“^) where ^ is the magnetic Fermi 
boundary associated with D = NIL. Also, since for any fixed a e Z 

N + a a 

lim - = Dp and lim — = 0 (4.29) 

L—>+oo L L—»+oo L 

one has 


lim XN+a ^ Qf 

L—>+oo 


and lim Xa = -Qf ■ (4.30) 

L ^+00 

These pieces of information are enough so as to obtain the large-L expansion of the particle-hole roots. One gets 


ttN+a - Qf + 


a-f^^\QF) 
^'.{QF I Qf) 


+ 


°(^) 


~ Qf "I" 


a-ff{-QF) 
^k{-QF I Qf) 


+ 




(4.31) 


Note that, above, the dependence of on the position of the particle-hole roots is kept implicit. Inserting the 
expansion of the roots up to 0(L“') one gets that 

I {7pJ1 ; {7hXf) - ^ - £{Z{A | g^) - i) + ^[i - ^(d, \ Q^) - ^(d, -Q^ \ Q^)] + o(i) . (4.32) 

Straightforward algebra based on these expansions leads to the claim. ■ 


Conclusion 

This paper develops tools allowing one to prove the existence and form of the large-L asymptotic expansion 
of the counting function associated with the XXZ spin-1/2 chain. The method is robust in that it does not rely on 
details of the model such as the sign or the magnitude of the L“ norm of the derivative of the bare phase. As such 
it allows one to step out of the setting where the model’s Yang-Yang action is strictly convex or where the model is 
a small perturbation of a non-interacting model. For these reasons it appears plausible that the method of large-L 
analysis proposed in this paper is applicable to many other quantum integrable models, this independently of the 
value of their coupling constant, viz. the sign of the Lieb kernels. 

The method should also allow one to study the case of excited states above the ground state described by 
Bethe roots such that a fixed number thereof is complex valued. I plan to investigate this issue in a forthcoming 
publication. 
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A Auxiliary results 

Below, by / holomorphic on F with F not open, it is meant that / is a holomorphic function on some open 
neighbourhood of F. Also, given a segment I the box-open neighbourhood is as defined in (13.101) while 

stands for the canonically oriented boundary of Srj^eiO- 

Proposition A.l, Let I = [a',b] be a segment and X be a compact in C such that I c Int(A). Let fi be a sequence 
of holomorphic functions on X converging to a holomorphic function f on X. Assume that f is such that sgn(/j|^ .^j) 
is constant. Then, there exists q,€ > 0 such that 

f ■ -S 2 T^pe{I) f{S 2 ri, 2 e{I)) is a biholomorpMsm. (A.l) 

Furthermore, there exists Lq such that for any L> Lq one has the inclusions 

Ml) C f{SMl)) C fL{S2r,,2eiI)) (A.2) 


and 


fl : S2,,2M)nf[\fiSr,M)))^fi'Sr,Ai)) 


(A.3) 


is a biholomorphism. Also, there exists a constant Urj^Af) depending on f such that, for any L > Lq, it holds 

inf inf \MA) - z| > uMf) > 0 ■ (A-4) 

AedS2r,.2AI) zef{Sr,,AI)) 


Finally, for any z e f{SrjM))’ tt holds 

^ IIA-ZIHx, (A.5) 

dS2,j.2dI) 

for some constant C[/] > 0 only depending on f, e, q and the segment 1. 


Proof — 

I first establish the statement relative to /. It is enough to consider the case where > 0. Then, by 

continuity, there exists Ty, e > 0 such that ^ f'W 0 for any A € Hence, / is a local 

biholomorphism on Syq^yM)- Thus it remains to show that there exists q,rf>q>0 such that i^ injective. 

If not, then there exists two sequences 

^ , a-3€< %{zf) <b + 3e such that 3(zl”^) ^ 0 and /(zf) - /(z^”^) • (A.6) 

I Z K K i Z 

By compactness of one may build converging subsequences zi^"^ —> Xa e [a - 3e\b + 3e\. Yet since 

f\[a-^e-b+^£\ > 0> / is bijective on \a - 3e',b + 3e] and thus f{x\) - f{x 2 ) implies that one has xi = X 2 = x. 
There exists an open neighbourhood Ux of x such that fu^ is a biholomorphism. For n large enough, one has that 
z^f"^ € Ux, a - 1,2. Thence, by injectivity of / on Ux, Zj^"^ - z^^"^ for n large enough, a contradiction. Thus / is a 
biholomorphism on Syj^M^) for some q > 0 and, a fortiori (lA.ll) is established. 

It remains to establish the statements relative to M Since / is injective on ii follows that 

f{dS 2 ^. 2 e{l)) Fl f{S,jM)) - ®- Thus, by compactness of these sets, one has 

2m;;, 6[/] ^ inf inf |/(z) -/(5)| > 0 . (A.7) 

zeSn,Al) sedS2„.2AI )' ' 
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Furthermore, the invertibility of / on S2rj,i£(J) ^Iso ensures that 


%,,[/] - inf(|/'(i)| : ^ £ 52,,2.(/)) > 0 . (A.8) 

Since /l —> / on X, one has that there exists Lq such that 

II/l - ^ min{M;,,f[^]/Ci, 4,;,[/]} provided that L > Lq . (A.9) 

Here, C\ is such that 


ll/'ll.~(5,,.(/)) ^ Cl • ll/'llL»(.53„3d/)) 


(A. 10) 


This bound implies (IA.4I) and ensures that 2i;, > 0 so that fi is a local biholomorphism on 

Sm, 2 e{l)- Furthermore, in virtue of Rouche’s theorem. 


for any z £ f{S„e{I)) s fiis) - z 


(A. 11) 


has exactly one simple zero in kS 2 ;,, 2 e(f)^ and is thus injective on this set. Therefore, it follows that 

A : S2,MI)^fE'{fi‘S,M))) ^ fiS.An) (A. 12) 

is a biholomorphism. The above also implies that c fL(<S2r/,2f;(^))- Also, since / is strictly increasing on 

[a - e; Z? + f] it is clear that for L large enough /([a - e',b + e\) D /l(/). From there the statement about inclusion 
follows. 

Finally, the bounds (IA.7I) and (IA.9I ) ensure that the integral representation (IA.5I) for A holds. The same repre¬ 
sentation holds with / /l so that 


fz\z)-r\z) 



dSiriMI) 


/i(2)(/(d)-z)-/'(d)(/L(T)-z) 

(AW - z)if{A) - z) 


lin 


The difference is then bounded thanks to (IA.4I) . (IA.7I) . (IA.9I) and the bound (lA.lOl) . 


(A. 13) 
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